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ABSTRACT 


This  report  summarizes  our  research  on  the  interaction  of  infrared  radiation, 
especially  from  C02  lasers,  with  the  cornea.  The  research  reported  here  was 
carried  out  between  July  1,  1977  and  June  30,  1985.  The  report  discusses  laser  sys¬ 
tems,  animals,  endothelial  temperature  histories,  endothelial  damage  thresholds, 
stromal  damage,  wound  repair,  epithelial  damage  from  single-  and  multiple- 
pulse  exposures,  and  a  new  thermal  damage  model  based  on  a  phase  transition. 
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SUMMARY 


This  report  summarizes  our  research  on  the  interac¬ 
tion  of  infrared  radiation,  especially  from  C02  lasers, 
with  the  cornea.  The  research  reported  here  was  con¬ 
ducted  between  July  1,  1977  and  June  30,  1985.  The 
following  topics  are  addressed: 

•  Laser  systems, 

•  Animals, 

•  Endothelial  temperature  histories, 

•  Endothelial  damage  thresholds, 

•  Stromal  damage, 

•  Wound  repair, 

•  Epithelial  damage  from  single-  and  multiple-pulse 
exposures,  and 

•  Thermal  damage  model. 

Sections  1.0  and  2.0  contain  general  descriptions  of 
the  laser  systems,  methods  of  beam  characterization, 
and  animal  care  and  handling  methods  used  through¬ 
out  the  reporting  period.  More  specialized  experimen¬ 
tal  details  are  presented  as  necessary  in  the  appropriate 
sections. 

In  Section  3.0  we  report  calculations  of  the  endothelial 
temperature  caused  by  exposing  the  cornea  to  C02  la¬ 
ser  radiation.  The  calculations  assume  that  the  laser  is 
operated  with  a  Gaussian  beam  intensity  profile  and  that 
the  radiation  is  absorbed  according  to  the  Beer-Lambert 
law.  The  efficacy  of  the  theoretical  model  is  tested  by 
comparing  the  calculations  with  measurements  on  ex¬ 
cised  corneas  exposed  to  C02  laser  radiation.  Both  a 
liquid  crystal  technique  and  a  thermocouple  were  used 
to  measure  the  temporal  and  spatial  variations  of  en¬ 
dothelial  temperature.  The  experimental  and  theoreti¬ 
cal  values  are  in  reasonable  accord,  which  gives  added 
confidence  in  the  theoretical  model,  especially  as  it  could 
be  applied  to  other  laser  systems  or  exposure  conditions. 

Section  4.0  reports  our  determinations  of  endotheli¬ 
al  damage  thresholds  in  rabbit  corneas  for  exposure  to 
C02  laser  radiation.  Endothelial  damage  was  detect¬ 
ed  using  a  staining  technique.  Threshold  damage  is  char¬ 
acterized  by  distorted  cells  and  uneven  staining  of  the 
cell  borders,  while  at  exposure  times  of  up  to  90%  of 
the  threshold,  the  stained  endothelia  were  indistinguish¬ 
able  from  those  in  control  eyes.  Three  peak  irradiance 
levels  (24.5,  10.0,  and  2.6  W/cm2)  were  investigated 
and  the  exposure  duration  was  varied  until  threshold 
damage  occurred.  In  each  instance  the  duration  of  the 
threshold  exposure  for  endothelial  damage  (1.0,  5.2,  and 
240  s,  respectively)  was  about  10  or  more  times  the  ac¬ 
cepted  duration  for  threshold  epithelial  damage.  These 
results  suggest  that  protecting  against  epithelial  dam¬ 
age  also  will  ensure  against  endothelial  damage  for 


C02  laser  systems.  However,  the  calculated  peak  tem¬ 
perature  increases  of  the  endothelium  at  its  threshold 
are  similar  to  those  of  the  epithelium  at  its  damage 
threshold.  We  show  that,  either  for  a  larger  beam  di¬ 
ameter  with  C02  radiation  or  for  other  lasers  where 
the  absorption  length  of  their  radiation  is  comparable 
to  the  corneal  thickness,  the  endothelial  temperature  rise 
can  approach  that  of  the  epithelium  for  the  same  ex¬ 
posure.  Thus,  under  these  conditions,  the  endothelial 
damage  threshold  could  be  similar  to  the  epithelial  dam¬ 
age  threshold. 

In  Section  5.0  we  note  that  exposures  well  above  the 
threshold  for  epithelial  damage,  but  below  that  for  endo¬ 
thelial  damage,  produce  bowl-shaped  stromal  wounds. 
In  such  wounds,  light  and  electron  microscopy  and  slit- 
lamp  photographs  all  show  a  sharp  demarcation  be¬ 
tween  the  damaged  and  undamaged  regions  48  h  after 
exposure.  The  micrographs  also  show  that  the  damaged 
region  is  acellular.  Calculations  of  the  expected  tem¬ 
perature  increases  combined  with  analyses  of  slit-lamp 
photographs  show  that  the  wound  boundary  corre¬ 
sponds  to  a  surface  of  equal  peak  temperature  increase. 
Comparisons  with  epithelial  and  endothelial  damage 
conditions  suggest  that  stromal,  endothelial,  and  epi¬ 
thelial  cells  have  essentially  the  same  thermal  damage 
mechanism. 

In  Section  6.0  we  discuss  the  healing  process  follow¬ 
ing  exposures  that  produce  stromal  wounds  extending 
about  halfway  through  the  cornea.  The  lesions  were 
documented  by  slit-lamp  photography  and  examined 
by  light  and  electron  microscopy.  Examinations  were 
carried  out  at  intervals  from  1  h  postexposure  up  to  1 1 
months  postexposure.  The  exposures  did  not  damage 
the  epithelial  basement  membrane.  The  cornea  is  cov¬ 
ered  by  a  mature  and  well-adhering  epithelium  by  1 
week.  The  stromal  wounds,  which  are  nearly  acellular 
at  48  h,  are  repopulated  with  active  fibroblasts  by  2 
weeks.  The  wounds  are  nearly  healed  by  8  weeks.  No 
inflammatory  response  or  vascularization  was  observed 
at  any  time  during  the  healing  period. 

In  Section  7.0  we  present  data  for  damage  thresholds 
of  the  corneal  epithelium  resulting  from  single-  and 
multiple-pulse  exposures.  To  our  knowledge,  this  is  the 
first  study  of  multiple-pulse  damage  to  the  cornea.  We 
confirm  that  damage  from  single-pulse  exposures  hav¬ 
ing  durations  between  10  3  and  10  s  can  be  correlat¬ 
ed  empirically  by  a  modified  critical  temperature  law 
in  which  the  critical  temperature  has  a  weak  dependence 
on  exposure  duration.  As  an  additional  test  of  the  modi¬ 
fied  critical  temperature  model,  we  also  examined  the 
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effect  of  beam  diameter  on  thresholds  for  single-pulse 
exposures.  The  results  supported  the  modified  critical 
temperature  model;  i.e.,  for  a  fixed  exposure  duration, 
the  threshold  irradiance  was  that  required  to  produce 
a  fixed  peak  temperature  increase,  independent  of  beam 
diameter.  The  damage  data  for  multiple  pulses  cover 
exposures  from  pulse  trains  of  up  to  999  pulses,  hav¬ 
ing  pulse  repetition  frequencies  between  1  and  100  Hz 
and  individual  pulse  durations  between  10  3  and  0.5 
s.  These  data  can  be  correlated  empirically  either  by  an 
approximate  critical  temperature  model,  by  a  model  in 
which  the  threshold  energy  density  varies  approximately 
as  the  }A  power  of  the  number  of  pulses,  or  by  a  mod¬ 
el  in  which  the  threshold  energy  density  varies  approx¬ 
imately  as  the  0.7  power  of  the  total  duration  of  the 
pulse  train. 


In  Section  8.0  we  introduce  a  new  physical  model  to 
explain  the  weak  dependence  of  the  critical  tempera¬ 
ture  on  exposure  duration  that  occurs  in  the  empirical 
modified  critical  temperature  model.  We  were  led  to 
this  new  model  by  energy  considerations  suggesting  that 
threshold  damage  can  be  associated  with  an  endother¬ 
mic  phase  transition  occurring  at  a  fixed  temperature. 
In  the  model,  we  associate  damage  with  a  fixed  amount 
of  energy  being  supplied  to  the  transition.  The  new  mod¬ 
el  can  correlate  all  of  the  single-pulse  threshold  data. 
It  predicts  a  transition  temperature  of  68°C  and,  if  we 
assume  that  approximately  5%  of  the  epithelium  un¬ 
dergoes  the  phase  change,  a  latent  heat  of  about  80 
cal/g.  The  analysis  suggests  that  the  weak  dependence 
of  critical  temperature  on  exposure  duration  in  the  em¬ 
pirical  model  actually  is  caused  by  heat  conduction  away 
from  the  damaged  region  during  exposure. 


FOREWORD 


In  conducting  the  research  described  in  this  report, 
the  investigators  adhered  to  the  Guide  for  the  Care  and 
Use  of  Laboratory  Animals,  prepared  by  the  Commit¬ 
tee  on  Care  and  Use  of  Laboratory  Animals  of  the  In¬ 
stitute  of  Laboratory  Animal  Resources,  National  Re¬ 


search  Council  (DHEW  Publication  No.  (N1H)  73-23, 
Revised  Edition,  1978). 

Citations  of  commercial  organizations  and  trade 
names  herein  do  not  constitute  an  official  Department 
of  the  Army  endorsement  or  approval  of  the  products 
or  services  of  these  organizations. 
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1.0  LASER  SYSTEMS 


Two  C02  lasers  were  used  during  the  contract  peri¬ 
od.  One  was  built  by  us  at  The  Johns  Hopkins  Univer¬ 
sity  Applied  Physics  Laboratory  (APL)  in  the  first  year 
as  partial  fulfillment  of  our  contract. 1  It  was  used  in 
all  endothelial  and  stromal  experiments  and  also  in  many 
of  the  epithelial  damage  threshold  experiments  for  both 
single-  and  multiple-pulse  exposures.  However,  the 
APL-built  laser  had  insufficient  power  output  to  pro¬ 
duce  damage  at  exposure  durations  less  than  10  ms.  We 
therefore  purchased  an  Apollo  Model  517  for  the  ex¬ 
periments.  Both  lasers  were  operated  in  the  TEMoq 
mode  for  which  the  beam  irradiance  follows  a  Gauss¬ 
ian  profile.  Laser  power  levels  were  measured  with  a 
Scientech  power  meter.  The  beam  irradiance  profile  was 
measured  at  each  experimental  session  with  either  a 
160-^m  circular  aperture  or  slits  having  widths  of  42 
or  16  pm.2  The  1/e  radius  of  the  beam  is  the  off-axis 
distance  at  which  the  irradiance  in  the  Gaussian  pro¬ 
file  is  1/e  (36.8%)  of  its  value  at  the  center.  For  a  Gauss¬ 
ian  beam  profile,  the  peak  irradiance,  70,  is  related  to 
the  total  power,  P,  in  the  beam  by  /0  =  P/A,  where 


A  is  the  area  within  the  1/e  radius.  The  nominal  1/e 
beam  diameter  of  the  Apollo  laser  was  4  mm,  while  that 
of  the  laser  built  at  APL  was  2.5  mm.  In  experiments 
requiring  different  beam  diameters,  a  zinc  selenide  lens 
with  a  12.5-cm  focal  length  was  inserted  into  the  opti¬ 
cal  system  and  the  beam  profile  was  determined  at  the 
position  where  the  cornea  would  be  located.  Damage 
thresholds  are  stated  as  the  peak  irradiance  or  peak  ener¬ 
gy  density  for  the  exposure  duration  under  investigation. 

Two  shutter  systems  were  used  to  control  the  exposure 
durations.  A  Uniblitz  shutter  was  used  for  both  single- 
and  multiple-pulse  experiments  in  which  individual  pulse 
durations  were  9  ms  or  greater.  The  shutter  was  oper¬ 
ated  by  a  control  system  that  allowed  us  to  set  the  du¬ 
ration,  number,  and  repetition  rate  of  the  pulses.  For 
pulses  shorter  than  9  ms,  a  chopper  wheel  was  used  in 
conjunction  with  the  Uniblitz  shutter.  In  this  configu¬ 
ration  the  chopper  wheel  controlled  the  individual  pulse 
duration  and  repetition  rate,  while  the  Uniblitz  was  used 
as  an  overall  gate  to  control  the  number  of  pulses. 


2.0  ANIMALS 


New  Zealand  white  rabbits  weighing  5  to  7  lb  were 
used  for  the  experiments.  During  the  first  several  years 
of  our  contract  we  administered  sodium  pentobarbital 
(Nembutal®)  to  anesthetize  the  animals.  In  particular, 
Nembutal®  was  employed  for  all  endothelial  and  stromal 
damage  experiments  (see  Sections  4.0  and  5.0).  Subse¬ 
quently  we  discovered  that  a  40:60  mixture  of  ketamine 
hydrochloride  (100  mg/ml)  and  xylazine  (20  mg/ml)  was 
a  more  controllable  general  anesthetic  for  rabbits  and 
also  was  better  for  stopping  eye  motion — an  especially 
useful  feature  with  longer  exposures  The  mixture  was 

'C.  B.  Bargeron,  W.  R.  Green,  R.  A.  Farrell,  and  R.  L.  McCally, 
•‘Structural  Alterations  in  the  Cornea  from  Exposure  to  Infrared 
Radiation.”  Annual  Report  to  Letterman  Army  Institute  of  Re¬ 
search,  Presidio  of  San  Francisco,  Calif.  (Jul  1978). 

:R.  L.  McCally,  “Measurement  of  Gaussian  Beam  Parameters," 
Appl.  Opt.  23,  2227  (1984). 


used  for  both  the  wound  healing  and  the  epithelial  dam¬ 
age  experiments  (Sections  6.0  and  7.0).  In  addition  to 
either  of  the  general  anesthesias,  proparacaine  hydroch¬ 
loride  ophthalmic  solution  (Alcaine®)  was  applied  top¬ 
ically  to  the  cornea.  The  anesthetized  animals  were 
placed  in  a  conventional  holder  and  were  positioned  with 
the  aid  of  a  He-Ne  alignment  laser  whose  beam  was 
colinear  with  the  C02  laser  output.  The  incident  radi¬ 
ation  was  aligned  perpendicular  to  the  surface  of  the 
cornea.  It  was  sometimes  necessary  to  use  a  speculum 
to  hold  the  eye  open,  particularly  for  the  longer  ex¬ 
posures.  In  such  instances  the  speculum  was  inserted 
approximately  1  min  before  exposure.  In  all  experi¬ 
ments,  the  cornea  was  irrigated  about  20  s  before  ex¬ 
posure  with  a  small  amount  of  physiological  saline  at 
room  temperature.  At  -10  s  before  exposure,  excess 
saline  was  blotted  by  holding  an  absorbent  tissue  against 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPUEO  PHYSICS  LABORATORY 

LAUREL.  MARYLAND 


the  limbus  just  below  the  area  to  be  exposed.  This  pro-  that  the  cornea’s  surface  temperature  returned  to  its 
cess  assured  a  reproducible  “tear”  layer.  Because  of  the  steady  state  value.  Following  exposure,  the  eye  was 

small  amount  of  liquid  involved,  the  10-s  interval  be-  blinked  manually  to  reform  its  natural  tear  film, 

tween  blotting  and  exposure  was  sufficient  to  ensure 


3.0  ENDOTHELIAL  TEMPERATURE  HISTORIES 


An  understanding  of  the  conditions  required  for  min¬ 
imal  endothelial  damage  as  discussed  in  Section  4.0  re¬ 
quires  knowledge  of  the  temperature-time  history  of  the 
endothelium  following  laser  exposure.  Consequently  we 
developed  methods  to  measure  and  to  calculate  the  spa¬ 
tial  distribution  of  temperature  at  the  position  of  the 
endothelium  as  a  function  of  time  following  exposure 
of  the  cornea’s  epithelial  surface  to  CO,  laser  radia¬ 
tion.  The  calculations  assume  that  the  cornea’s  thermal 
properties  are  homogeneous  throughout  its  thickness 
and,  for  convenience,  are  approximated  by  those  of  wa¬ 
ter.'4  The  measurements  use  a  well-characterized  C02 
laser  and  serve  to  test  the  theoretical  model.  Once  veri¬ 
fied  in  this  manner,  the  model  can  be  used  to  predict 
the  effects  of  other  laser  systems.5'6 

The  calculations  are  based  on  a  Green’s  function  so¬ 
lution  to  the  heat  flow  equation  for  a  Gaussian  inci¬ 
dent  beam.  The  beam  is  assumed  to  impinge  on  a  semi¬ 
infinite  slab  and  to  have  a  constant  peak  irradiance  for 
the  duration  of  the  exposure.  The  solution  T(r,z,t), 
where  r  is  the  radial  distance  from  the  beam  axis,  z  is 
the  depth  from  the  front  surface,  and  t  is  the  time  af¬ 
ter  the  initiation  of  exposure,  is  given  in  terms  of  a  def- 

'  V  S.  Brownell  and  B.  E.  Stuck,  “Ocular  and  Skin  Hazards  from 
C02  Laser  Radiation,”  in  Proc.  9th  Army  Science  Conf.,  pp. 
123-138  (Jun  1974). 

4D.  E.  Egbert  and  E.  F.  Maher,  Cornea I  Damage  Thresholds  for 
Infrared  laser  Exposure:  Experimental  Data,  Model  Predictions, 
and  Safety  Standards,  U.S.A.F.  School  of  Aerospace  Medicine, 
Brooks  AFB,  SAM-TR-77-29  (1977). 

5C.  B.  Bargeron,  R.  L.  McCally,  and  R.  A,  Farrell,  "Calculated 
and  Measured  Endothelial  Temperature  Histories  of  Excised  Rab¬ 
bit  Corneas  Exposed  to  Infrared  Radiation,”  Exp.  Eye  Res.  32, 
241-250  (1981). 

AC.  B.  Bargeron,  R.  A.  Farrell,  W.  R.  Green,  and  R.  L.  McCally, 
"Corneal  Damage  from  Exposure  to  Infrared  Radiation:  Rabbit 
Endothelial  Damage  Thresholds,"  Health  Phys.  40,  855-862 
(1981). 


inite  integral.17  In  Ref.  1  we  outlined  the  derivation 
of  this  integral,  both  for  the  case  of  total  surface  ab¬ 
sorption  and  for  Beer’s  law  absorption,  and  gave  a  list¬ 
ing  of  the  computer  program  (available  on  request)  that 
we  wrote  for  its  evaluation. 

This  theoretical  model  neglects  heat  loss  via  convec¬ 
tion,  radiation,  and  evaporation  from  the  warm  anterior 
surface.  We  showed  that  these  losses  are  insignificant 
for  the  exposure  conditions  encountered  in  the  en¬ 
dothelial  damage  threshold  experiments.5  (Since  these 
exposures  are  more  severe  than  those  in  the  stromal 
damage  experiments  discussed  in  Sections  5.0  and  6.0, 
it  follows  that  such  approximations  are  valid  there  as 
well.) 

The  Green’s  function  solution  does  not  account  for 
the  latent  heat  associated  with  phase  transitions  such 
as  tear-film  boiling,  which  could  occur  in  exposures 
where  the  anterior  surface  reaches  100°C,  or  protein 
denaturation,  which  occurs  at  about  60°C.4  Although 
methods  are  available  for  including  phase  transitions 
in  heat  flow  calculations  (see  Section  8.0), 7,8  we  elect¬ 
ed,  for  the  purposes  of  this  section,  to  use  the  Green’s 
function  approach  to  obtain  an  upper  bound  on  the  en¬ 
dothelial  temperature.  One  objective  of  the  experimental 
measurements  was  to  test  how  well  this  bound  approx¬ 
imates  the  actual  temperature. 

Experimental  temperatures  were  obtained  with  excised 
corneas  using  either  a  thermocouple  or  a  liquid  crystal 
technique. 5  The  thermocouple  measures  temperatures 
over  a  wide  range,  but  only  at  a  single  position  in  space. 
On  the  other  hand,  the  liquid  crystal  technique  has  a 
limited  practical  operating  range  and  is  less  precise,  but 

H.  S.  Carslaw  and  J.  C.  Jaeger,  Conduction  of  Heat  in  Solids, 

Oxford  University  Press,  Oxford  (1959). 

8S.  G.  Bankoff,  “Heat  Conduction  or  Diffusion  with  Change  of 

Phase,"  Adv.  Chem.  Eng.  5,  51-122  (1964). 
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measures  temperature  over  the  entire  endotheliaJ  surface. 

In  the  liquid  crystal  technique,  crystals  were  painted 
onto  thin  (25-^m)  black  mylar  disks  that  were  preformed 
to  match  the  corneal  curvature.  The  disks  were  attached 
at  their  edges  to  excised  corneas  using  cyanoacrylate  ce¬ 
ment.  The  corneas  were  mounted  in  an  appropriate 
holder  and  motion  pictures  were  made  of  the  color 
changes  in  the  liquid  crystals  as  the  cornea  responded 
to  the  laser  exposure.  The  films  were  subsequently  evalu¬ 
ated  as  described  in  Ref.  5  to  obtain  the  temperature- 
time  histories  at  various  radial  distances  from  the  laser 
beam  axis. 

Representative  data  obtained  using  the  liquid  crystal 
method  are  shown  in  Figs.  1  and  2  for  two  different 
exposures,  each  near  the  epithelial  damage  threshold 
for  the  corresponding  experimental  irradiance  level.  The 
error  bars  primarily  reflect  the  difficulty  in  estimating 
the  color  variation  of  the  crystals.  Figure  1  shows  the 
spatial  variation  of  the  temperature  increase  and  indi¬ 
cates  that  the  temperature  rise  falls  off  sharply  at  the 
1/e  radius  of  the  beam  (r1/e  =  1.0  mm).  Figure  2 
shows  the  temperature  increase  on  the  beam  axis  and 
compares  it  to  the  cal  elated  value. 
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Figure  1  The  measured  endothelial  temperature  in¬ 
crease  as  a  function  of  time  at  three  radial  positions 
from  the  beam  axis.  The  liquid  crystal  technique  was 
used  for  these  measurements.  The  time,  f  =  0,  corre¬ 
sponds  to  the  initiation  of  exposure.  The  peak  irradiance 
was  11.0  W/cm2,  the  exposure  duration  was  0.50  s,  and 
the  1/e  radius  of  the  beam  was  1.0  mm.  The  circles  rep¬ 
resent  the  data  on  the  beam  axis  (r  =  0);  the  squares 
represent  the  data  at  r  =  0.5  mm;  and  the  triangles  rep¬ 
resent  the  data  at  r  =  1.0  mm.  The  error  bars  indicate 
an  uncertainty  of  ±0.5°C. 


Figure  2  The  circles  are  liquid  crystal  measurements 
of  the  endothelial  temperature  increase  on  the  beam  axis 
as  a  function  of  time.  The  curve  is  the  calculated  value. 
The  peak  irradiance  was  25.1  W/cm2,  the  duration  was 
0.104  s,  and  the  1/e  radius  of  the  beam  was  1.0  mm.  For 
the  calculation,  the  thickness  was  taken  as  0.500  mm 
to  account  for  the  corneal  thickness  measured  at  zero 
transcorneal  pressure  and  the  thickness  of  the  plastic 
film. 

The  thermocouple  temperature  measurements  were 
made  using  a  0.025-mm-diameter  chromel-alumel  ther¬ 
mocouple.  The  small  thermal  inertia  associated  with  this 
thermocouple  gave  it  a  rapid  response  time,  and  the  con¬ 
ductive  heat  loss  from  the  junction  via  the  small  di¬ 
ameter  wires  was  negligible  for  our  experimental  condi¬ 
tions.  (Details  of  constructing  and  mounting  the  ther¬ 
mocouple  are  given  in  Ref.  5.)  The  junction  was  posi¬ 
tioned  carefully  at  the  endothelial  surface  of  the  mounted 
cornea  with  the  aid  of  a  slit-lamp  microscope.  Accurate 
positioning  of  the  probe  on  the  center  of  the  laser  beam 
was  accomplished  by  moving  the  cornea  holder/thermo¬ 
couple  assembly  with  an  x-y  micropositioner  until  we 
obtained  the  maximum  temperature  rise  for  a  very  low 
exposure  (i.e.,  25  W/cm2  for  0.05  s,  which  is  well  be¬ 
low  the  epithelial  damage  threshold).  This  method  en¬ 
sured  a  positioning  accuracy  to  within  about  0. 1  mm 
for  the  2-mm-l /e-diameter  beam. 

Representative  temperature-time  histories  of  the  en¬ 
dothelium  that  were  measured  using  the  thermocouple 
are  shown  in  Fig.  3,  wherein  peak  irradiance  was  24 
W/cm2.  The  exposure  durations  in  Figs.  3a  and  3b 
were,  respectively,  nearly  equal  to  and  10  times  the 
epithelial  damage  threshold.  The  latter  exposure  was 
chosen  because  it  is  near  the  threshold  for  endothelial 
damage6  (cf.  Section  4.0).  The  curves  in  Fig.  3  are  the 
values  calculated  from  the  model.  The  different  tem¬ 
perature  increases  in  Fig.  3a  and  in  Fig.  2,  both  cor- 
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Figure  3  (a)  The  circles  are  measurements  of  the  en¬ 
dothelial  temperature  increase  as  a  function  of  time 
made  with  the  0.025-mm  thermocouple  positioned  on 
the  beam  axis.  The  curve  is  the  value  calculated  for  the 
same  position.  The  peak  irradiance  was  24.0  W/cm2, 
the  duration  was  0.104  s,  and  the  1/e  radius  of  the  beam 
was  1.0  mm.  The  measured  in  situ  comeal  thickness  was 
0.38  mm.  (b)  The  measured  endothelial  temperature  in¬ 
crease  as  in  (a)  with  the  exception  that  the  exposure  dur¬ 
ation  was  increased  to  1.04  s.  This  exposure  corre¬ 
sponds  to  the  threshold  for  endothelial  damage  at  this 
power  density.6 

responding  to  similar  exposures,  can  be  attributed  to 
the  different  corneal  thicknesses  in  the  two  experiments. 
If  the  thickness  difference  is  taken  into  account,  the  two 
measuring  techniques  are  in  close  agreement. 

Figures  2  and  3  show  that  the  measured  temperatures 
are  consistently  somewhat  lower  than  the  temperatures 
calculated  from  the  model.  Nevertheless,  the  agreement 
is  quite  reasonable.  Several  factors  associated  with  the 
temperature  measurements  (e.g.,  imperfect  thermal  con¬ 
tact  between  film  or  thermocouple  and  the  cornea,  fi¬ 
nite  response  times,  heat  loss  in  the  thermocouple  leads, 
etc.),  which  were  assumed  to  be  small,  all  act  to  reduce 
the  measured  temperature.  Also,  as  noted  above  in  the 


discussion  of  the  theoretical  model,  the  neglect  of  heat 
losses  from  the  anterior  surface  via  evaporation,  con¬ 
vection,  and  radiation,  as  well  as  the  latent  heats  as¬ 
sociated  with  phase  transitions,  serves  to  make  the  calcu¬ 
lations  represent  an  upper  bound  to  the  actual  endotheli¬ 
al  temperature.  Thus  the  good  agreement  achieved — 
even  in  the  exposure  of  Fig.  3a,  where  all  of  these  ef¬ 
fects  would  have  been  most  severe — justifies  the  approx¬ 
imations  of  the  model. 

Another  approximation  is  the  neglect  of  possible  con¬ 
vective  heat  transfer  at  the  endothelial  surface.  Indeed, 
we  did  observe  convection  in  the  sample  cell  while  mak¬ 
ing  the  liquid  crystal  temperature  measurements.9  The 
Green’s  function  solution  assumes  that  the  conduction 
equation  holds  across  the  endothelial  boundary  and  that 
the  cornea  and  bathing  media  have  the  same  thermal 
properties.  Convection  currents  would  change  the  semi¬ 
infinite  slab  problem  into  a  finite  thickness  slab  prob¬ 
lem  with  a  convective  boundary  condition  at  the  en¬ 
dothelial  surface.  10  This  boundary  condition  relates 
the  temperature  difference  across  the  boundary  to  its 
normal  derivative  there.  Conventionally,  one  assumes 
a  linear  proportionality  (which  introduces  an  unknown 
constant),  and  uses  an  eigenfunction  expansion  to  cal¬ 
culate  the  temperature  increases. 10  We  did  develop  the 
expansion  appropriate  for  both  the  uniform  beam  prob¬ 
lem  and  the  Gaussian  beam  problem,9  and,  as  expect¬ 
ed,  the  solutions  showed  a  more  rapid  reduction  in 
endothelial  temperature  after  attainment  of  the  peak 
value  than  for  the  more  slowly  cooling  conducting  case. 
However,  no  a  priori  procedure  is  available  for  deter¬ 
mining  the  unknown  proportionality  constant.  Although 
the  constant  could  be  “determined”  a  posteriori  so  as 
to  give  a  best  fit  to  the  experimental  data,  we  deem  such 
an  approach  arbitrary,  especially  in  view  of  the  other 
approximations  discussed  above.  Consequently,  we 
elected  to  ignore  possible  convection  currents  and  use 
the  semi-infinite  slab  approximation  for  our  calcula¬ 
tions,  particularly  since  this  approximation  is  in  such 
good  agreement  with  the  experiments. 


"C  B.  Bargeron.  R.  A.  Farrell.  \V.  R.  Green,  and  R  L.  McCally. 
"Structural  Alterations  in  the  Cornea  from  Exposure  to  Infrared 
Radiation."  Annual  Report  to  l  etlerman  Army  Institute  of  Re¬ 
search.  Presidio  of  San  Francisco,  Calif.  (Sep  1979). 

"H.  V  Olcer.  "On  the  Theory  of  Conductive  Heat  Transfer  in 
Finite  Regions  with  Boundary  Conditions  of  the  Second  Kind,” 
Ini  J  Hem  Mass  Transfer  8.  529-556  (1965). 
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4.0  ENDOTHELIAL  DAMAGE  THRESHOLDS 


This  research  stemmed  from  a  report  that  there  was 
some  evidence  of  altered  endothelial  integrity  follow¬ 
ing  laser  exposures  at,  or  slightly  above,  the  threshold 
level  for  epithelial  damage. n  If  correct,  this  observa¬ 
tion  could  have  clinical  implications  because,  in  primates 
and  humans,  endothelial  damage  is  not  believed  to  be 
repaired  by  cell  division,  but  rather  by  the  undamaged 
cells  enlarging  and  sliding  to  cover  the  wound. 11  Thus 
the  effect  of  repeated  exposures  could  be  cumulative 
and  eventually  the  endothelium’s  ability  to  maintain  nor¬ 
mal  corneal  hydration  (and  therefore  transparency) 
might  be  impaired. 

Endothelial  damage  thresholds  were  determined  for 
peak  irradiances  of  24.5,  10.0,  and  3.6  W/cm.6  The 
C02  laser  was  operated  as  usual  in  the  TEMqo  mode 
with  a  1/e  beam  diameter  of  2.0  mm.  Exposure  times 
were  increased  from  values  near  the  (then)  known  epi¬ 
thelial  damage  thresholds  until  we  found  the  minimum 
time  required  to  produce  endothelial  damage. 

To  detect  endothelial  damage,  the  rabbits  were  sac¬ 
rificed  with  an  intravenous  injection  of  Nembutal®  ap¬ 
proximately  2  h  postexposure.  The  globes  were  enu¬ 
cleated  and  the  corneas  were  excised  at  the  limbus.  Dam¬ 
age  was  detected  by  wet  staining  the  endothelium  with 
alizarin  red  S  and  trypan  blue  dyes.”  In  some  instances 
indocyanine  green  dye  was  used  in  place  of  trypan  blue 
with  equivalent  results. 14  In  this  staining  method  aliza¬ 
rin  red  S  dye  clearly  delineates  the  borders  of  both 
healthy  and  damaged  cells,  while  trypan  blue  dye  (or 
indocyanine  green  dye)  penetrates  those  cells  whose 
membranes  have  become  leaky. 

The  damage  thresholds  for  the  three  irradiance  lev¬ 
els  are  shown  in  Fig.  4. 6  Results  were  obtained  on 
several  rabbits  at  each  exposure  level.  The  average  thick¬ 
ness  of  the  corneas  was  about  0.4  mm.  The  same  ex¬ 
posure  duration  was  found  for  threshold  damage  when 
the  corneas  were  examined  either  24  or  48  h  after  the 
10-W/cnT  exposure;  thus,  there  are  no  apparent  long¬ 
term  effects  beyond  those  evident  at  2  h. 

"E.  S.  Beatrice  and  B.  E.  Stuck,  “Ocular  Effects  of  Laser  Radia¬ 
tion:  Cornea  and  Anterior  Chamber,’’  AGARD  Lecture  Series 
Vol.  79,  Laser  Hazards  and  Safety  in  the  Military  Environment, 
p.  5  (1975). 

I:D.  L.  Van  Horn  and  R.  A.  Hyndirck,  “Endothelial  Wound  Re¬ 
pair  in  Primate  Cornea,"  Exp.  Eye  Res.  21,  113-124  (1975). 
mD.  J.  Spence  and  G.  A.  Peyman,  “A  New  Technique  for  the  Vi¬ 
tal  Staining  of  the  Corneal  Endothelium,"  Invest.  Ophthalmol. 
15,  1000-1002  (1976). 

I4J.  K.  McEnerney  and  G.  A.  Peyman,  “Indocyanine  Green— A 
New  Vital  Stain  for  Use  Before  Penetrating  Keratoplasty,"  Arch. 
Ophthalmol.  %,  1445-1447  (1978). 
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Figure  4  Peak  irradiance  versus  endothelial  damage 
threshold  exposure  time.  The  1/e  radius  of  the  beam  was 
1.0  mm  for  all  exposures.  The  circles  indicate  the  mini¬ 
mum  time  for  which  we  observed  damage  and  the  par¬ 
tial  error  bar  extends  to  the  longest  exposure  for  which 
no  damage  was  observed.  Sixteen  rabbits  were  used  to 
determine  the  threshold  at  3.6  and  10  W/cm2;  five  rab¬ 
bits  were  used  at  the  24.5-W/cm2  condition. 


Two  hours  after  the  cornea  is  exposed  at  the  en¬ 
dothelial  damage  threshold,  slit-lamp  examination  rev¬ 
eals  severe  corneal  edema  localized  at  the  lesion  and 
strong  stromal  light  scattering  to  a  depth  of  about  V* 
to  Vs  the  corneal  thickness.  The  severity  of  the  epithelial 
lesion  and  damage  to  the  anterior  stroma  at  these  ex¬ 
posures  depend  strongly  on  the  peak  irradiance.  At  24.5 
W/cm:  (1.04  s),  the  epithelium  is  destroyed  and  there 
is  significant  stromal  cratering,  whereas  at  3.6  W/cm2 
(240  s),  the  epithelium  is  damaged,  but  not  destroyed, 
and  no  stromal  dimpling  or  cratering  is  evident.  These 
observations  are  not  surprising  because  calculations 
show  that  the  thermal  insult  to  the  anterior  cornea  is 
much  less  severe  for  the  3.6-W/cm:  (240-s)  exposure 
than  it  is  for  the  higher  irradiance,  but  shorter  dura¬ 
tion,  exposures. 

The  appearance  of  the  endothelium  for  various  ex¬ 
posure  conditions  is  illustrated  in  Figs.  5,  6,  and  7.  Fig¬ 
ure  5  shows  the  endothelium  of  a  cornea  2  h  after 
exposure  to  10-W/cm:  radiation  for  about  90S  of  the 
endothelial  threshold  damage  time  (i.c. ,  for  4.56  s).  The 
cells  appear  like  those  in  control  corneas  and  no  dam¬ 
age  is  evident.  This  result  also  was  found  in  corneas 
examined  24  and  48  h  after  similar  exposures.  The  un¬ 
even  focusing  across  the  photograph  is  typical  for  these 
wet  preparations,  which  arc  often  not  completely  flat. 
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Figure  5  A  photograph  of  endothelial  cells  near  the 
beam  axis.  The  peak  irradiance  was  10  W/cm2,  the  du¬ 
ration  was  14.56  s  (corresponding  to  about  0.9,  the  dam¬ 
age  threshold  time),  and  the  1/e  radius  of  the  beam  was 
1.0  mm.  The  endothelial  cells  appear  normal. 
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Figure  6  A  low  magnification  photograph  of  the  en¬ 
dothelial  surface  of  a  cornea  exposed  to  10-W/cm2  radi¬ 
ation  for  5.24  s.  The  1/e  radius  of  the  beam  was  1.0  mm. 
This  exposure  was  near  the  endothelial  damage  thresh¬ 
old.  The  outer  circular  ring  defines  a  portion  of  the  dam¬ 
aged  epithelial  surface  as  seen  through  the  preparation. 
The  inner  circular  region  is  the  endothelial  thermal  le¬ 
sion.  The  dark  indocyanine  green  staining  between 
points  A  and  B  is  associated  with  a  folding  artifact  and 
not  with  thermal  damage.  The  arrow  points  to  the  re¬ 
gion  shown  in  Fig.  7. 
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Figure  7  A  higher  magnification  view  of  the  region  in¬ 
dicated  by  the  arrow  in  Fig.  6. 


Figure  6  is  a  low-magnification  view  of  the  cornea’s 
posterior  surface  2  h  after  irradiation  at  10  W/cm2  for 
5.24  s,  which  is  at  (actually  slightly  above)  the  endotheli¬ 
al  damage  threshold.  The  extent  of  anterior  corneal 
damage  also  can  be  seen  through  such  preparations.  At 
threshold  exposures  for  endothelial  damage,  this  view' 
of  the  anterior  damage  usually  shows  two  distinct  con¬ 
centric  areas,  the  inner  area  probably  being  associated 
with  the  aforementioned  stromal  cratering.  The  endo¬ 
thelial  lesion  is  represented  by  the  circular  region  at  the 
center  of  Fig.  6.  The  center  of  this  lesion  is  not  coinci¬ 
dent  with  that  of  the  anterior  surface  damage  area  (prob¬ 
ably  as  a  result  of  the  preparation).  Threshold  endotheli¬ 
al  damage  is  characterized  by  distorted  cells  and  uneven 
staining  of  the  cell  borders.  In  particular,  much  of  each 
damaged  cell  border  is  more  heavily  stained  with  alizarin 
red  S  than  are  normal  cell  borders;  yet  other  parts  of 
the  damaged  cell  borders  are  often  completely  unstained. 
Figure  7  highlights  these  features.  The  location  of  the 
area  depicted  in  Fig.  7  is  indicated  by  the  arrow  in  Fig. 
6.  Another  common  feature  of  these  lesions  is  that  the 
boundary  between  the  damaged  area  and  the  region  con¬ 
taining  apparently  normal  cells  is  rather  sharp  (com¬ 
pare  Figs.  6  and  7). 

At  exposures  that  are  somew  hat  above  threshold,  en¬ 
dothelial  cells  at  the  center  of  the  lesion  are  obliterated 
and  cells  immediately  surrounding  the  bare  region  ap¬ 
pear  to  have  been  pushed  away  from  the  center.  Further 
outward,  the  cells  look  similar  to  those  in  a  threshold 
lesion,  and  as  with  a  threshold  lesion,  there  is  a  sharp 
boundary  with  the  undamaged  area. 
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We  found  that  alizarin  red  S  stain  gives  the  most  use¬ 
ful  damage  information  for  thermal  lesions  of  the  en¬ 
dothelium.  Trypan  blue  or  indocyanine  green  staining 
is  sporadic  and  sometimes  absent.5  When  present,  the 
staining  is  frequently  associated  with  preparation  ar¬ 
tifacts  such  as  folds  (cf.  Fig.  6). 

The  endothelial  temperature  increases  on  the  beam 
axis  were  calculated  for  the  threshold  conditions  of  Fig. 
4  and  are  shown  as  functions  of  time  in  Fig.  8.  At 
threshold,  the  endothelium  experiences  peak  tempera¬ 
ture  elevations  of  30  to  50°C.  These  values  are  similar 
to  those  experienced  by  the  epithelium  at  its  damage 
threshold3,411  (see  also  Section  7.0).  Thus  the  endothe¬ 
lium’s  sensitivity  to  thermal  insult  is  similar  to  that  of 
the  epithelium. 

The  data  in  Fig.  4  show  that  for  a  C02  laser  with 
a  2-mm-l/e  beam  diameter,  the  exposure  durations  re¬ 
quired  to  cause  threshold  endothelial  damage  are  about 
10  times  those  required  to  produce  threshold  epithelial 
damage  at  the  same  irradiances.  This  result  would  not 
necessarily  hold  for  other  lasers  or  even  for  other  ex¬ 
posure  conditions  with  the  C02  laser.  Figure  9  illus¬ 
trates  this  point  for  the  C02  laser  (wavelength,  10.6 
fim),  the  holmium  laser  (wavelength,  2.06  jim),  and  the 
erbium  laser  (wavelength,  1.54  /cm).5’6  Differences  in 
the  results  for  the  three  lasers  are  caused  by  the  marked 
variations  in  the  corresponding  absorption  length  of 
their  radiation  in  the  cornea  (assumed  to  be  the  same 
as  in  water).  (The  absorption  length  is  the  inverse  of 
the  absorption  coefficient  in  the  exponential  expression 
of  the  Beer-Lambert  law.  It  is  the  distance  into  the  medi¬ 
um  at  which  63.2%  of  the  incident  radiation  has  been 
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Figure  8  Calculated  endothelial  temperature  histories 
for  exposures  of:  a)  24.5  W/cm2,  1.04  s;  b)  10  W/cm2, 
5.24  s;  and  c)  3.6  W/cm2,  240  s.  At  240  s,  curve  c) 
reaches  about  32'C  as  indicated  by  curve  c')  and  cor¬ 
responding  time  scale  along  the  top  of  the  figure.  A  cor¬ 
neal  thickness  of  0.40  mm  and  a  1/e  beam  radius  of  1 .0 
mm  were  was  used  in  the  computation.  Zero  time  cor¬ 
responds  to  the  beginning  of  the  exposure. 


absorbed.)  The  absorption  lengths  for  holmium  and  er¬ 
bium  lasers  are,  respectively,  256  and  840  pm,15  or 
about  one-half  and  twice  the  corneal  thickness;  the  value 
for  the  C02  laser  is  only  10.6  /xm,  or  about  one-fourth 
the  thickness  of  the  epithelium.4 


Exposure  time  (s) 


Figure  9  (a)  Predicted  temperature  increase  as  a  func¬ 
tion  of  exposure  time  for  a  TEM,*,  mode  C02  laser  beam 
of  peak  irradiance  1  W/cm2.  Note  that  since  the  temper¬ 
ature  increase  is  directly  proportional  to  the  peak  irradi¬ 
ance,  the  increases  for  other  values  of  peak  irradiance 
can  be  found  by  multiplication.  The  solid  curves  repre¬ 
sent  temperature  increases  at  a  depth  of  10  /xm  (anterior 
epithelium);  the  dotted  curves,  at  a  depth  of  400  pm  (en¬ 
dothelium  for  rabbit).  The  value  of  d,  the  1/e  diameter  of 
the  beam,  is  in  millimeters.  The  absorption  coefficient 
is  950  cm  “ 1  for  10.6-/xm  radiation,  (b)  Same  as  (a)  except 
the  conditions  are  for  a  holium  laser  (wavelength,  2.06 
pm;  absorption  coefficient,  39  pm  "’).  (c)  Same  as  (a)  ex¬ 
cept  the  conditions  are  for  an  erbium  laser  (wavelength, 
1.54  pm;  absorption  coefficient,  11.9  pm  "’). 

’’K.  F.  Palmer  and  D.  Williams,  ‘‘Optical  Properties  of  Water  in 
the  Near  Infrared,"  J.  Opt.  Soc.  Am.  64,  1107-1110  (1974). 
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Figure  9  shows  that  the  ratio  of  the  exposure  time 
needed  to  cause  a  given  endothelial  temperature  increase 
to  the  exposure  time  needed  to  cause  the  same  increase 
in  epithelial  temperature  will  decrease  with  increasing 
beam  diameter  and/or  with  increasing  absorption  length. 
For  these  curves  the  ratio  is  as  small  as  1.25  for  the 
erbium  laser  and  as  large  as  10  for  a  2-mm-diameter 
CO,  laser  beam.  The  results  in  Fig.  9  also  show  that 
the  ratio  of  the  epithelial  temperature  increase  to  the 


endothelial  temperature  increase  at  a  fixed  exposure  time 
decreases  with  increasing  exposure  time  but  is  indepen¬ 
dent  of  the  peak  irradiance.  Thus,  for  a  given  epitheli¬ 
al  temperature  increase,  the  concomitant  endothelial 
temperature  increase  will  be  larger  with  a  lower  peak 
irradiance  beam.  Therefore,  the  characteristics  of  the 
infrared  laser  beam  must  be  specified  carefully  when 
discussing  potential  hazards  to  endothelial  cells. 


5.0  STROMAL  DAMAGE 


During  the  investigation  of  endothelial  damage  we 
noted  that  corneas  exposed  at  greater  than  the  epithelial 
damage  threshold,  but  less  than  the  endothelial  dam¬ 
age  threshold,  develop  bowl-shaped  lesions  in  the  stro¬ 
ma  that  are  nearly  devoid  of  cells  48  h  postexpo¬ 
sure.161  At  that  time,  such  corneas  have  recovered 
their  normal  thickness  and  have  regrown  a  smooth 
epithelium.  Figure  10,  a  montage  of  light  micrographs, 
clearly  shows  this  acellular  region.  The  border  between 
the  damaged  and  undamaged  region  is  sharply  delineat¬ 
ed.  The  abruptness  of  the  demarcation  suggested  to  us 
that  keratocyte  damage,  like  that  of  endothelial  and 
epithelial  cells,  depends  sensitively  on  the  thermal  insult. 

To  investigate  this  dependence,  we  exposed  several 
corneas  so  that  the  damage  extended  about  halfway 
through  the  stroma.  Two  exposure  levels  were  inves¬ 
tigated,  one  at  9.7  W/cm2  (peak  irradiance)  for  2.5  s 
and  the  other  at  26  W/cm2  for  0.4  s.  The  1/e  beam 
diameters  were  1 .88  mm  for  the  first  exposure  level  and 
2.44  mm  for  the  second.  The  durations  of  both  ex¬ 
posures  are  about  four  times  longer  than  the  duration 
required  to  produce  minimal  epithelial  damage3'4'1  IS 


l6R.  L.  McCally,  C.  B.  Bargeron,  W.  R.  Green,  and  R.  A.  Farrell. 
"Stromal  Damage  in  Rabbit  Corneas  Exposed  to  CO,  Laser 
Radiation,"  Exp.  Eye  Res.  37,  543-550  (1983). 

1  R.  A.  Farrell,  R.  L.  McCally,  C.  B.  Bargeron,  and  W.  R.  Green, 
“Alterations  in  the  Cornea  from  Exposure  to  Infrared  Radiation,” 
Annual  Report  to  Letterman  Army  Institute  of  Research,  Presi¬ 
dio  of  San  Francisco,  Calif.  (Sep  1982). 

'"N.  A.  Peppers,  A.  Vassiliadis,  L.  G.  Dedrick,  H.  Chang,  R.  R. 
Peabody,  H.  Rose,  and  H.  C.  Zweng,  "Cornea  Damage 
Thresholds  for  CO.  Laser  Radiation,”  Appl.  Opt.  8,  377-381 
()%9). 


(see  Section  7.0)  and  both  exposures  produced  similar 
damage  to  the  cornea.  We  took  slit-lamp  photographs 
within  1  h  after  the  exposures  and  again  at  48  h  after 
which  the  rabbits  were  sacrificed  and  the  exposed  eyes 
were  prepared  for  light  and  electron  microscopy  (see 
Ref.  17  for  details).  (A  few  animals  were  also  sacrificed 
at  1  h.) 

Figure  1 1  shows  the  appearance  of  the  cornea  at  1 
and  48  h  after  an  exposure  to  26  W/cm2  for  0.4  s.  At 
1  h,  the  slit  lamp  reveals  a  circular  hazy  lesion  whose 
periphery  is  surrounded  by  a  raised  ring  of  epithelium. 
As  noted  above,  by  48  h  the  cornea  has  regained  its 
normal  thickness,  its  front  surface  is  covered  by  newly 
grown  epithelial  cells,  and  it  is  completely  smooth.  The 
circular  hazy  lesion  remains.  The  bowl-shaped  lesion 
is  clearly  visible  in  the  narrow-slit  photograph. 

Figure  10,  which  was  referred  to  above,  is  a  com¬ 
posite  light  micrograph  of  a  cornea  exposed  like  the  one 
in  Fig.  1 1 .  The  sections  were  cut  as  near  as  possible  to 
a  plane  passing  through  the  center  of  the  circular  dam¬ 
age  area.  The  nearly  cell-free  lesion  extends  about  half¬ 
way  through  the  stroma  at  its  deepest  point.  Transmis¬ 
sion  electron  micrographs  also  were  made  of  this  le¬ 
sion.  Figure  12a  shows  that  the  basement  membrane 
is  intact.  At  this  stage  of  healing,  however,  the  basal 
epithelium  is  flattened  and  irregular.  The  anterior  stro¬ 
ma  shows  a  greatly  reduced  number  of  keratocytes. 
Those  that  remain  appear  necrotic  with  dark,  dense, 
intracytoplasmic  deposits,  small  vacoules,  and  a  loss  of 
normal  intracytoplasmic  organelles.  The  sharpness  of 
the  demarcation  between  damaged  and  undamaged  stro¬ 
ma  is  clearly  revealed  at  the  higher  magnification  of  Fig. 


5§y 


Figure  12  (a)  A  low-magnification  electron  micrograph 
of  the  anterior  stroma  and  epithelium  of  the  cornea 
shown  in  Fig.  10.  The  basal  epithelium  is  flattened  and 
irregular,  but  the  basement  membrane  is  intact.  The  ker- 
atocytes  that  remain  appear  necrotic  with  dark,  dense, 
intracytoplasmic  deposits,  small  vacoules,  and  a  loss 
of  normal  intracytoplasmic  organelles,  (b)  A  higher  mag¬ 
nification  electron  micrograph  of  the  same  cornea  show¬ 
ing  the  sharpness  of  the  lesion  edge.  A  necrotic  kera- 
tocyte  lies  only  8  anterior  to  a  completely  normal 
keratocyte. 


12b,  which  shows  a  necrotic  keratocyte  only  8  ^m  an¬ 
terior  to  a  completely  normal  keratocyte. 

We  used  high-magnification  slit-lamp  photographs, 
taken  just  before  sacrificing  the  animal  for  histology, 
to  measure  the  shape  of  the  lesion.  The  slit  was  made 
very  narrow  for  these  photographs  and  the  lamp  was 
arranged  so  that  the  light  was  incident  normal  to  the 
cornea’s  surface.  The  narrow,  illuminated  region  was 
positioned  so  that  it  passed  directly  through  the  center 


of  the  circular  lesion,  allowing  us  to  obtain  a  cross  sec¬ 
tion  of  the  brightly  scattering  bowl-shaped  lesion  at  its 
deepest  point.  The  position  of  the  lesion  border  was 
measured  from  the  photographs  and  the  depths  were 
corrected  for  refraction  at  the  anterior  corneal  surface. 
Figure  13  shows  the  border  positions  on  two  corneas 
with  different  exposures.  The  curves  are  calculated 
isotherms.  For  both  exposures  the  lesion  border  close¬ 
ly  follows  the  isotherms. 
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Figure  13  Measured  position  of  the  border  of  a  bowl¬ 
shaped  stromal  lesion.  The  points  are  measured  from 
a  highly  magnified  slit-lamp  photograph  taken  48  h  after 
an  exposure  of  (a)  9.7  W/cm2  for  2.5  s  with  a  1/e  beam 
radius  of  0.94  mm  and  (b)  26  W/cm2  for  0.4  s  with  a  1/e 
beam  radius  of  1.21  mm.  The  curves  are  calculated 
isotherms. 


The  calculations  also  reveal  that  the  temperature  his¬ 
tory  is  similar  at  different  positions  on  a  given  peak  tem¬ 
perature  isotherm.  This  is  illustrated  in  Fig.  14,  which 
shows  the  temperature  histories  for  the  9.7-W/cm2- 
2.5-s  exposures  at  a  point  on  the  beam  axis  and  a  depth 
of  243  j*m,  and  at  a  point  having  a  radial  position  695 
jim  from  the  beam  axis  and  a  depth  of  50  jim.  These 
two  positions  are  both  on  the  45°C-peak-temperature 
isotherm.  As  expected,  the  temperature  rise  lags  slight¬ 
ly  for  the  deeper  position,  but  the  shape  of  the  two 
curves  is  nearly  identical.  Similar  results  were  obtained 
at  other  positions  for  this  exposure  and  also  for  the 
26-W/cm:-0.4-s  exposure. 

These  experiments  show  that  the  border  of  the 
stromal  lesion  produced  by  CO;  laser  irradiation  is 
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Figure  14  Temperature  histories  at  two  positions  on 
the  isothermal  surface  having  a  45°C  temperature  rise 
following  a  9.7-W/cm2-2.5-s  exposure.  The  1/e  radius  of 
the  beam  was  taken  as  0.94  mm.  The  solid  curve  is  on 
the  beam  axis  at  a  depth  of  243  jim.  The  dashed  curve 
is  at  a  point  695  jtm  from  the  beam  axis  at  a  depth  of 
50  jim. 


sharply  delineated  (as  viewed  in  both  the  slit  lamp  and 
in  the  histologic  sections)  and  that  it  corresponds  to  a 
surface  of  equal  peak  temperature  rise.  The  conditions 
on  this  border  represent  the  threshold  for  stromal  cel¬ 
lular  damage,  which  is  at  a  slightly  higher  temperature 
for  shorter  exposure  durations  than  it  is  for  longer  ex¬ 
posure  durations.16 

Similar  peak  temperatures  and  dependencies  on  ex¬ 
posure  duration  have  been  found  for  epithelial  and  en¬ 
dothelial  damage3,4'617'18  (also  see  Sections  4.0  and 
7.0).  Thus  we  conclude  that  all  corneal  cells  have  a  simi¬ 
lar  thermal  damage  mechanism. 


6.0  WOUND  REPAIR 


As  discussed  above,  that  corneas  exposed  above  the 
threshold  for  epithelial  damage,  but  below  that  for  endo¬ 
thelial  damage,  develop  a  sharply  defined  bowl-shaped 
lesion,  with  borders  corresponding  to  a  surface  of  equal 
peak  temperature  increase  within  the  stroma.16'1  The 
peak  temperature  value  depends  on  exposure  duration, 
being  slightly  lower  for  longer  exposures.  Since  similar 


peak  temperatures  and  dependencies  have  been  found 
for  threshold  epithelial  and  endothelial  dam¬ 
age,'46'1  18  we  concluded  that  all  corneal  cells  have  a 
similar  thermal  damage  mechanism.  Since  the  basement 
membrane  did  not  seem  to  be  damaged  in  these  experi¬ 
ments,  '6''  at  least  up  to  48  h  postexposure  (see  Sec¬ 
tion  5.0),  the  time  course  of  healing  for  this  type  of 
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wound  was  of  considerable  interest.  Consequently,  we 
extended  our  work  on  such  exposures  to  document  the 
healing  process. w 

All  exposures  were  made  on  the  central  cornea  us¬ 
ing  a  1/e  beam  diameter  of  2.8  mm.  The  irradiance  was 
fixed  at  25  W/cm:  and  the  exposure  duration  was  0.4 
s.  As  noted  in  Section  5.0,  this  exposure  level  produces 
a  wound  that  penetrates  about  halfway  through  the  cor¬ 
nea  at  its  deepest  point.  In  the  first  series  of  experiments, 
14  rabbits  were  exposed.  Slit-lamp  photographs  were 
taken  at  1,  24,  and  48  h  and  at  1,2,  4,  and  8  weeks 
postexposure.  Two  rabbits  were  sacrificed  at  each  in¬ 
terval  and  the  exposed  corneas  were  prepared  for  histo¬ 
logical  examination.  Representative  slit-lamp  photo¬ 
graphs  and  histological  results  are  shown  in  Figs.  15 
through  25.  In  a  second  set  of  experiments  we  exposed 
six  more  rabbits.  Three  of  the  rabbits  were  used  for  slit- 
lamp  observation  at  3,  4,  and  5  days  and  one  eye  was 
prepared  for  histological  e.\,.mination  at  each  of  these 
times.  The  remaining  three  rabbits  were  used  for  slit- 
lamp  observations  up  to  1 1  months  postexposure. 

Fig.  15a  shows  the  appearance  of  the  cornea  at  1  h. 
The  center  of  the  wound  is  highly  scattering,  the  cen¬ 
tral  cornea  is  thickened  and  rough,  and  the  epithelium 
is  raised  around  the  periphery  of  the  lesion.  There  is 
a  sharp  transition  from  the  damaged  area  to  the  clear, 
evidently  undamaged,  cornea.  (I  ig.  15a  is  similar  to  Fig. 
11.)  F  igure  16  shows  typical  light  micrographs  of  the 
wound  area.  In  F  ig.  16a,  the  damaged  area  extends 
about  halfway  through  the  stroma  and  also  indicates 
the  irregular  roughening  of  the  central  ephithelium.  The 
endothelium  and  IX’scemet’s  membrane  appear  normal. 
The  transition  from  normal  to  burned  epithelium  is 
shown  in  Fig.  16c.  The  raised  area,  which  is  evident 
in  the  slit-lamp  photograph,  .s  clearly  indicated.  Inside 
this  area  the  epithelium  is  thinned  and  darkly  stained 
(see  also  Fig.  16b).  In  the  electron  microscope  the 
epithelium  shows  marked  disorganization,  including  nu¬ 
clear  chromatin  clumping,  loss  of  plasma  membranes, 
and  loss  ot  cytoplasmic  organelles  (F  ig.  19).  The  base¬ 
ment  membrane  looks  undamaged.  F  igure  20  shows  a 
necrotic  keratocyte  in  the  anterior  stroma  and  abnor¬ 
mally  fuzzy  or  smudged  collagen  fibrils. 

The  appearance  of  the  cornea  at  24  h  is  shown  in 
Figs.  15b  and  15c.  In  some  cases  at  24  h,  the  burned 
epithelium  is  completely  opaque  and  has  not  yet  peeled 

'’R.  A.  harrell,  C  B  Bargcron.  R.  I  Met  ally,  and  W  R.  Green, 
"Structural  Alterations  in  the  Cornea  from  Exposure  to  Infrared 
Radiation,"  Annual  Report  to  l  etterman  Army  Institute  of  Re¬ 
search,  Presidio  of  San  I  rancisco,  Calif.  (Jan  1 98S ) 

:"R.  I  Viet  ally.  C.  B  Bargcron.  VV  R  Green,  and  R  A.  I  arrell. 
"Beam  Diameter  Dependence  and  Healing  Processes  in  CO,  1  a 
ser  Damaged  C  orneas,"  Ahs  Suppl  .  Imoi  Ophthalmol  l  isa 
a I  Sa.  25.  .128  (1984). 


off  (cf.  Fig.  15b),  whereas  in  other  cases,  the  burned 
epithelium  has  disappeared,  revealing  a  brightly  scat¬ 
tering  wound  in  the  stroma  (cf.  Fig.  15c,  except  that 
here,  two  small  ‘‘flakes”  of  opaque  material  remain  in 
the  central  area).  Figure  17  shows  that  the  wound  is  al¬ 
ready  covered  by  an  immature  epithelial  layer  that  is 
only  about  30  gm  thick.  The  immature  epithelium  con¬ 
sists  of  basal  cells,  which  have  not  yet  attained  their 
usual  columnar  form,  covered  by  a  single  layer  of  super¬ 
ficial  cells.  These  structures  are  seen  in  more  detail  in 
ihe  electron  micrograph  of  Fig.  21,  which  also  shows 
some  debris  of  keratocytes  in  the  anterior  stroma. 

The  acellular  nature  of  the  stromal  wound  at  48  h 
(discussed  extensively  in  Refs.  16  and  17  in  Section  5.0) 
is  confirmed  by  the  light  micrograph  shown  in  Fig.  17b. 
This  micrograph  also  shows  that  the  epithelium  now 
has  several  layers  of  wing  cells  and  that  the  columnar 
cells  are  better  formed.  The  electron  micrograph  in  Fig. 
22  shows  the  new  epithelium  in  more  detail  and  reveals, 
as  noted  before,  that  the  anterior  stroma  is  devoid  of 
cells.  Only  the  debris  of  a  few  keratocytes  remains. 

The  appearance  of  the  lesion  at  1  and  2  weeks  is 
shown  in  Figs.  15d  and  15e.  At  1  week,  the  lesion  looks 
very  much  as  it  does  at  48  h.  By  2  weeks,  however,  the 
anterior  surface  of  the  cornea  is  slightly  flattened  as  a 
result  of  stromal  thinning  (cf.  Fig.  15e)  and  the  lesion 
is  less  visible  in  the  wide-slit  photograph  of  Fig.  15e. 
Histology  reveals  that  the  epithelium  is  completely  nor¬ 
mal  at  1  week  (Fig.  17c)  and  that  the  wound  area  is 
still  acellular.  The  few  dark  staining  areas  are  appar¬ 
ently  debris  of  the  necrotic  keratocytes  shown  by  the 
electron  microscope  (Fig.  23).  By  2  weeks,  however,  the 
stroma  is  repopulated  by  keratocytes  (Figs.  24  and  25) 
and  the  stromal  collagen  has  regained  its  normal  ap¬ 
pearance  (cf.  Fig.  25). 

The  stromal  thinning  and  anterior  flattening  are  still 
evident  at  4  and  8  weeks  (cf.  Fig.  15f  and  15g).  At  4 
weeks  there  is  still  some  scattering,  but  it  is  confined 
to  the  anterior-most  portion  of  the  stroma.  The  lesion 
remains  faintly  visible  in  the  wide-slit  photograph  of 
F  ig.  !5f  and  is  clearly  visible  in  narrow-slit  photographs. 
At  8  weeks  the  anterior  scattering  has  diminished  but 
is  visible  in  the  narrow -slit  photograph  of  Fig.  15g.  The 
hazy  lesion  is  almost  invisible  in  the  wide-slit  photo¬ 
graph  of  Fig.  15g.  The  light  micrographs  show  that  the 
keratocytes  that  have  repopulated  the  lesion  area  con¬ 
tinue  to  stain  darker  than  normal  (Fig.  18). 

No  involvement  of  vascular  components  (e.g.,  mono¬ 
cytes  or  polymorphonuclear  leukocytes)  or  any  vascular¬ 
ization  was  observed  during  the  experiments.  Also,  the 
endothelium  and  Desccmet’s  membrane  remained  nor¬ 
mal  throughout  the  8-vveck  period.  The  second  series 
of  experiments  at  days  3,  4,  and  5  also  failed  to  show 
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any  evidence  of  inflammatory  cells,  which  could  have 
been  missed  in  the  first  series  of  experiments.  No  re¬ 
markable  evidence  was  uncovered  as  the  second  series 


u 


Figure  15  Slit-lamp  photographs  of  laser-induced  le¬ 
sions  produced  by  a  24-W/cm2  peak  irradiance  and  a 
0.4-s  duration.  The  1/e  radius  of  the  beam  was  1.4  mm, 
(a)  At  1  h  postexposure,  a  raised  ring  on  the  periphery 
of  the  lesion  is  typical,  (b)  At  24  h  postexposure,  the  le¬ 
sion  is  covered  by  the  completely  opaque,  burned 
epithelium  that  has  not  yet  sloughed  off.  About  one-half 
of  the  rabbits  have  this  appearance  at  24  h.  (c)  Also  at 
24  h  postexposure,  the  opaque  covering  has  nearly 
sloughed  off.  There  are  two  small  “flakes"  remaining  in 
the  center,  (d)  1  week  postexposure,  (e)  At  2  weeks 
postexposure,  the  stroma  has  thinned  slightly  in  the  cen¬ 
ter  of  the  lesion,  causing  a  flattening  of  the  anterior  cor¬ 
neal  surface,  (f)  At  4  weeks  postexposure,  stromal 
thinning  and  anterior  flattening  are  now  quite  evident. 
The  stromal  scattering  is  significantly  reduced  and  is 
confined  to  the  anterior-most  region  of  the  stroma,  (g) 
At  8  weeks  postexposure,  the  central  area  of  the  lesion 
remains  thin  and  the  ha^y  lesion  is  only  faintly  visible. 
Photographs  (a),  (e),  (f)  and  (g)  are  all  of  the  same  rabbit. 


of  rabbits  was  observed  at  8  months  and  again  at  1 1 
months  (the  slit-lamp  photographs  of  these  rabbits  up 
to  8  weeks  were  similar  to  those  of  the  first  set  of  rab¬ 
bits).  By  8  months,  however,  it  was  impossible  to  see 
the  wound  area  using  wide-slit  illumination.  In  the 
narrow-slit  view  a  very  slight  increase  in  scattering  re¬ 
mained  in  the  anterior-most  region  of  the  stroma.  The 
stroma  remained  thin  in  this  region,  as  had  first  been 
noted  at  2  weeks.  The  general  appearance  of  the  cor¬ 
neas  of  these  three  rabbits  was  virtually  the  same  at  1 1 
months  as  it  was  at  8  months. 
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Figure  16  Light  micrographs  of  laser-induced  lesions  after  receiving  exposures  like  those  for  the  slit- 
lamp  photographs  (cf.  Fig.  15).  The  1-jim-thick  Epon  sections  were  stained  with  hematoxylin  and  eosin. 
At  1  hour,  (a)  the  center  of  the  lesion  extends  about  halfway  through  the  stroma.  The  central  epithelium 
is  highly  disrupted  and  roughened,  (b)  The  epithelium  near  the  center  is  coagulated  but  smooth,  (c)  The 
edge  of  the  lesion  showing  the  transition  between  undamaged  and  damaged  epithelium.  The  raised  and 
disrupted  edge  that  was  evident  in  the  slit  lamp  (cf.  Fig.  15a)  is  shown  clearly.  The  damaged  epithelium 
stains  very  darkly  and  is  significantly  thinner,  although  all  of  the  cell  layers  are  evident. 
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Figure  17  Same  as  Fig.  16.  (a)  The  central  area  of  the  lesion  at  24  h.  The  lesion  is  covered  by  a  new 
epithelial  layer;  however,  its  basal  cells  have  not  attained  their  columnar  shape  and  there  is  only  a  sin¬ 
gle  layer  of  flat  wing  cells,  (b)  The  central  area  of  the  lesion  at  48  h.  The  basal  cells  are  closer  to  their 
usual  columnar  form  and  there  are  several  layers  of  wing  cells,  (c)  The  central  area  of  the  lesion  at  1 
week.  The  epithelium  is  normal  in  appearance.  A  few  dark  staining  areas  in  the  anterior  stroma  are  evi¬ 
dently  debris  of  damaged  keratocytes.  (d)  The  central  area  of  the  lesion  at  2  weeks.  The  lesion  is  repop¬ 
ulated  with  keratocytes,  which  appear  to  stain  more  darkly  than  usual. 
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Figure  18  Same  as  Figs.  16  and  17.  The  central  area  of  the  lesion  at  (a)  4  and  (b)  8  weeks.  The  kerato- 
cytes  continue  to  stain  darkly. 
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Figure  19  Electron  micrograph  of  the  epithelium  and  anterior  stroma  in  the  central  lesion  area  at  1  h 
(cf.  Figs.  15a  and  16).  The  epithelium  shows  marked  disorganization,  including  nuclear  chromatin  clumping, 
loss  of  plasma  membranes,  and  loss  of  cytoplasmic  organelles  The  basement  membrane  appears  to 
be  undamaged. 
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Figure  20  Electron  micrograph  of  the  anterior  stroma  in  the  central  lesion  area  at  1  h  showing  a  necrotic 
keratocyte  (cf.  Figs.  15a,  16.  and  19).  Note  also  the  abnormal  fuzzy  or  smudged  appearance  of  the  colla¬ 
gen  fibrils. 


Figure  21  Electron  micrograph  of  the  epithelium  in  the  central  lesion  area  at  24  h  (cf.  Figs  15b.  15c. 
and  17a).  The  new  epithelium  consists  of  only  two  layers  of  cells  The  basal  cells  have  not  attained  their 
usual  columnar  form  and  there  is  only  a  single  layer  of  flat  wing  cells  This  micrograph  also  shows  the 
debris  of  a  few  keratocytes  in  the  anterior  stroma 
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Figure  22  Electron  micrograph  of  the  epithelium  and  anterior  stroma  in  the  central  lesion  area  at  48  h 
(cf.  Figs.  10.  lie,  and  1 7b).  There  are  several  layers  of  wing  cells  and  the  basal  cells  are  better  formed 
than  at  24  h  (cf.  Fig.  21).  The  anterior  stroma  is  devoid  of  cells;  only  the  debris  of  necrotic  keratocytes 
remains. 


Figure  23  Electron  micrograph  of  the  epithelium  and  anterior  stroma  in  the  central  lesion  area  at  1 
week  (cf  Figs.  I5d  and  1 7c).  The  basal  cells  are  normal  and  the  wound  area  remains  acellular  except 
for  some  debris. 
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EPITHELIAL  DAMAGE  FROM  SINGLE-  AND  MULTIPLE-PULSE 

EXPOSURES 
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Figure  24  Electron  micrograph  of  the  epithelium  and 
anterior  stroma  in  the  central  lesion  area  at  2  weeks  (cf. 
Figs.  16e  and  17d).  The  stroma  is  repopulated  with  ker- 
atocytes. 


Figure  25  Electron  micrograph  of  the  basement  mem¬ 
brane  and  anterior  stroma  in  the  central  lesion  area  at 
2  weeks  (cf.  Figs.  15e,  17d,  and  24).  The  collagen  fibrils 
have  regained  their  normal  appearance. 


This  report  has  thus  far  dealt  with  thermal  damage 
to  the  endothelium  and  stroma  from  CO;  laser  radia¬ 
tion;  however,  most  of  the  work  of  others  has  concen¬ 
trated  on  epithelial  damage.  This  is  natural,  since  99°xi 
of  the  radiation  from  a  CO;  laser  is  absorbed  within 
the  first  50  gm  of  the  tear  film  and  epithelium.  The  ener¬ 
gy  is  rapidly  converted  to  heat  that  initially  is  concen¬ 
trated  in  the  volume  of  absorption.  Many  aspects  of 
CO-  laser  damage  to  the  corneal  epithelium  have  been 
investigated  over  the  past  15  years.  Egbert  and  Maher4 


summarized  the  earlier  work  and  discussed  the  empiri¬ 
cal  critical  temperature  and  damage  integral  models  for 
correlating  threshold  epithelial  damage.  The  data  they 
summarized  came  from  many  sources  and  covered  a 
wide  range  of  exposure  parameters  for  single-pulse  ex¬ 
posures  (multiple-pulse  data  were  not  then  availble  in 
the  literature).  Although  the  data  were  somewhat  scat¬ 
tered,  Egbert  and  Maher  found  that  both  of  these 
models  could  correlate  the  damage  conditions. 
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An  early  formulation  of  the  damage  integral  was 
made  by  Henriques,21  who  suggested  that  biological 
thermal  damage  depends  on  the  integrated  thermal  his¬ 
tory,  weighted  by  an  Arrhenius  rate  factor.  He  proposed 
the  damage  integral  in  the  form 


ft  =  A  J  exp(  - B/T(t))dt  ,  (1) 


1  where  A  and  B  are  constants  and  T(i)  is  the  absolute 

temperature  of  the  tissue  as  a  function  of  time.  Thresh¬ 
old  damage  is  reached  when  ft  =  1 .  The  constants  A 
and  B  are  to  be  determined  from  measurements  of  dam¬ 
age  thresholds.  The  model  was  applied  to  corneal  tis¬ 
sue  by  Brownell  and  Stuck.’  It  worked  reasonably  well 
for  the  earlier  data;  however,  the  values  they  found  for 
A  and  B  (A  =  exp(  141.6)  and  B  =  48,400  K)  were  such 
that  their  physical  significance  remained  obscure. 

The  critical  peak  temperature  ( CPT)  model  was  in- 
‘  troduced  by  Egbert  and  Maher.  In  it,  the  threshold  of 

|  damage  was  assumed  to  occur  when  the  tissue  was  heat- 

i  ed  to  a  certain  critical  temperature.  Egbert  and  Maher 

!  found  it  necessary  to  make  the  CPT  a  weak  function 

of  the  exposure  duration.  Because  the  CPT  is  not  con¬ 
stant  it  would  perhaps  be  more  accurate  to  say  that  they 
proposed  a  modified  peak  temperature  model  of  dam¬ 
age.  Nevertheless,  they  found  that  the  equation 

!  CPT  =  79.6  r  "'"5:  (2) 


gave  a  good  description  of  the  existing  data.  In  Eq.  2, 
the  CPT  is  in  degrees  centigrade  and  r  is  the  exposure 
duration  in  seconds. 

Both  the  damage  integral  and  CPT  models  depend 
on  being  able  to  determine  the  temperature  history  at 
some  position  in  the  epithelium,  which  is  usually  ac¬ 
complished  computationally  by  solv  ing  the  heat  equa¬ 
tion  discussed  in  Section  3.0. 

As  noted  above,  previous  studies  have  concentrated 
on  damage  caused  by  a  single  pulse  of  radiation  from 
a  CO;  laser.  In  some  applications,  however,  lasers  are 
operated  in  a  mode  in  which  they  emit  a  train  of  pulses. 
We  have  made  the  first  measurements  of  epithelial  dam¬ 
age  thresholds  for  exposures  to  sequences  of  subthresh¬ 
old  pulses.  We  also  obtained  data  for  single-pulse  expo¬ 
sures  with  durations  ranging  between  10  1  and  10  s 
using  the  same  lasers,  the  same  beam  conditions,  and 


I  C  Henriques,  Jr..  "Studies  of  thermal  Injurs  V.  the  Pre¬ 
dictability  and  the  Significance  of  Thermally  Induced  Rate  Pro¬ 
cesses  I  eadmg  to  Irreversible  Epidermal  Injury."  -1A/-I  Arch. 
Pathol.  4,1,  4X9  ( 194“), 
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the  same  damage  criterion  to  help  us  understand  the 
multiple-pulse  results.  In  addition,  we  measured  the  in¬ 
fluence  of  beam  diameter  on  single-pulse  damage  thresh¬ 
olds.  By  properly  interpreting  the  measurements,  we 
have  provided  corroboration  of  the  modified  critical 
temperature  model  for  correlating  damage.  Moreover, 
certain  aspects  of  the  results  with  single-pulse  exposures 
suggest  a  new  physical  model  that  can  explain  the  slight 
variation  in  the  CPT  with  exposure  duration,  at  least 
for  single-pulse  exposures.  The  new  model,  which  asso¬ 
ciates  thermal  damage  with  the  occurrence  of  an  endo¬ 
thermic  phase  transition,  is  discussed  in  Section  8.0. 

We  use  the  standard  criterion  for  minimal  epithelial 
damage,  which  is  the  presence  of  a  superficial  grayish 
white  spot  on  the  cornea  that  develops  within  0.5  h  of 
the  exposure.  ’ 4  Corneas  were  assessed  for  damage  by 
examination  with  a  Nikon*  photo  slit  lamp  about  30 
min  after  the  exposure.  For  near-threshold  damage,  the 
faint  spot  is  best  located  with  a  wide  slit.  Occasionally 
it  was  helpful  to  stain  with  fluorescein  as  a  secondary 
method  of  damage  assessment.  After  washing  away 
residual  stain,  any  thermal  damage  appears  as  a  dis¬ 
tinctive  circular  spot  of  diffuse  fluorescence. 

In  practice,  we  find  that  the  damage  threshold  is  shar¬ 
ply  defined,  and  we  therefore  do  not  use  techniques  such 
as  probit  analysis  for  its  determination.  Rather,  we  ex¬ 
pose  the  tissue  above  threshold,  then  below.  The  gap 
above  and  below  is  narrowed  to  bracket  the  threshold 
condition.  The  final  “bracket  width"  is  such  that  the 
threshold  value  of  irradiancc  is  known  to  a  precision 
of  approximately  ±5^0. 


7.1  SING1.E-PIT.SE  EXPOSE  RES 

The  experimental  results  for  exposure  durations,  r, 
between  -  10  ‘  and  —10s  are  summarized  in  Table 
1 .  The  table  also  shows  the  peak  temperature  rise  (above 
the  ambient  corneal  temperature)  on  the  beam  axis  10 
ixm  below  the  anterior  corneal  surface.  The  tempera¬ 
ture  was  calculated  using  the  Green's  function  solution 
discussed  in  Section  3.0.  Although  these  temperature 
rises  are  nearly  constant  (especially  in  the  range  10  ' 
<  r  <  1,  where  they  average  40  ±  2°C),  the  slight 
dependence  on  exposure  duration  noted  by  Egbert  and 
Maher  is  obvious  when  the  entire  range  is  taken  into 
account.  The  dependence  of  the  peak  temperature  as¬ 
sociated  with  damage  on  exposure  duration  is  exam¬ 
ined  further  in  l  ig.  26.  In  accordance  with  Egbert  and 
Maher,  we  have  plotted  the  natural  log  of  the  temper¬ 
ature  (in  degrees  centigrade),  and  not  the  temperature 
increase.  A  value  of  35 °C  is  assumed  for  the  ambient 
temperature  of  the  cornea's  anterior  surface.  A  least 
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Table  1 

Single-pulse  exposures. 


h 

(W/cm2)* 

Pulse  Duration 
(s) 

ED,h 

(J/cm') 

d 

(cm) 

fcr 

624 

0.00096 

0.599 

0.102 

54.0 

177 

0.0039 

0.690 

0.180 

47.0 

86.2 

0.0092 

0.790 

0.252 

41.4 

64.3 

0.0150 

0.964 

0.178 

42.2 

51.1 

0.0185 

0.945 

0.180 

38.1 

57.2 

0.0191 

1.09 

0.254 

43.6 

38.8 

0.0309 

1.20 

0.260 

39.5 

20.4 

0.101 

2.07 

0.158 

39.5 

8.99 

0.498 

4.48 

0.182 

37.9 

9.39 

0.500 

4.69 

0.224 

41.0 

6.46 

0.977 

6.31 

0.244 

38.3 

2.89 

9.730 

28.1 

0.248 

34.1 

‘No  errors  are  shown;  however,  the  values  are  assumed  to  carry  er¬ 
rors  of  ±  10%.  These  reflect  both  the  bracket  spacing  in  determining 
the  threshold  (±5%)  and  the  assumed  accuracy  of  the  power  mea¬ 
surement  (±5%).  The  same  percentage  errors  would,  of  course,  be 
reflected  to  ED,h  and  A7"max  . 

*  Maximum  temperature  rise  calculated  on  the  beam  axis  10  jim  be¬ 
neath  the  cornea’s  anterior  surface. 


squares  fit  yields  the  (empirical)  modified  critical  tem¬ 
perature  model  in  the  form 


CPT  =  72  t'm:  "C  ,  (3) 

which  is  in  close  agreement  with  the  correlation  obtained 
by  Egbert  and  Maher  (cf.  Eq.  2). 

Theory  predicts  that  the  irradiance  required  to  pro¬ 
duce  a  given  temperature  rise  for  a  given  single-pulse 
exposure  should  depend  on  the  diameter  of  the  beam. 
Therefore,  if  a  modified  critical  temperature  model  such 
as  that  in  Eq.  3  were  correct,  epithelial  damage 
thresholds  should  depend  on  beam  diameter  in  the  same 
way.  We  tested  this  prediction  by  measuring  damage 
thresholds  at  a  variety  of  beam  diameters.  Three  ex¬ 
posure  durations  were  used  for  the  experiments:  0.1, 
1.0,  and  10  s.  The  experimental  results  are  shown  by 
the  data  points  in  Fig.  27,  where  the  measured  thresh¬ 
old  irradiance  is  plotted  as  a  function  of  beam  diameter. 
Each  of  the  three  data  sets  shows  that  the  irradiance 
required  for  threshold  damage  increases  as  the  beam 


.1 


Exposure  duration  (s) 


Figure  26  The  natural  log  of  the  calculated  tempera¬ 
ture  10  jim  beneath  the  cornea's  anterior  surface  as  a 
function  of  exposure  duration.  The  ambient  temperature 
of  the  cornea  was  assumed  to  be  35°C.  The  ±  10%  er¬ 
ror  bars  reflect  both  the  bracket  spacing  in  determin¬ 
ing  the  threshold  and  the  estimated  accuracy  of  the 
power  measurement.  The  line  is  a  least  squares  fit  to 
the  data  and  yields  the  empirical  correlation  given  by 
Eq.  3. 
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diameter  decreases.  This  result  is  attributable  to  radial 
heat  conduction  in  the  cornea,  which  becomes  more  im¬ 
portant  as  the  irradiated  spot  size  decreases.  If  we  were 
to  assume  that  the  modified  CPT  model  were  correct, 
then  we  could  calculate  the  peak  irradiance  needed  to 
cause  this  temperature  rise  as  a  function  of  beam  di¬ 
ameter.  The  curves  in  Fig.  27  are  the  results  of  such 
calculations.  For  the  0.1-s  exposures,  the  curve  is  for 
the  peak  irradiances  needed  to  cause  a  41  °C  peak  tem¬ 
perature  rise  at  a  depth  of  10  ^m  along  the  beam  axis. 
The  curves  through  the  1.0-  and  10-s  exposure  data  are 
for  the  irradiances  needed  to  produce  peak  tempera¬ 
ture  increases  of  39°  and  36°C,  respectively.  In  all  cases 
the  peak  temperature  increases  are  close  to  the  values 
in  Table  1  for  the  respective  exposure  durations  and  are 
therefore  in  close  agreement  with  the  correlation  given 
in  Eq.  3. 

The  relationship  in  Eq.  3  is,  of  course,  only  an  em¬ 
pirical  correlation  of  the  damage  data.  To  gain  physical 
insight  into  the  basis  for  the  slight  dependence  of  the 
CPT  on  exposure  duration,  it  is  helpful  to  examine  how 
the  energy  density  needed  to  cause  minimal  damage, 
ED,h,  depends  on  exposure  duration.  The  values  of 
EDlh  from  Table  1  are  displayed  graphically  in  Fig.  28 
as  a  function  of  exposure  duration.  The  values  of 
ED,h  and  their  associated  exposure  durations  are  ex¬ 
perimental  data  (peak  energy  densities  of  the  Gaussian 
profiles)  and  are  therefore  independent  of  any  thermal 
model  assumed  for  temperature  calculations.  Our  inter¬ 
pretation  of  this  plot,  however,  assumes  either  that  all 
the  data  are  for  the  same  diameter  laser  beam  or  that 
the  data  are  obtained  under  conditions  for  which  the 
irradiance  required  to  produce  damage  is  independent 
of  beam  diameter  (cf.  Fig.  27).  At  very  short  exposures, 
ED,h  appears  to  be  asymptotically  approaching  a  fixed 
value.  This  behavior  might  be  expected  if  the  occurrence 
of  an  endothermic  phase  transition  within  the  epithelium 
were  the  underlying  damage  mechanism.  In  that  case,  a 
thermally  isolated  epithelium  would  sustain  damage 
when  it  absorbed  enough  energy  both  to  raise  its  temper¬ 
ature  to  the  transition  temperature  and  then  to  supply 
the  latent  heat  needed  to  cause  the  phase  change  (e.g., 
protein  denaturation,  alterations  of  the  organization  of 
the  lipids  in  the  cell  membranes,  etc.).  For  a  thermally 
isolated  epithelium,  ED,h  would  be  constant  and  inde¬ 
pendent  of  the  exposure  duration.  Very  short  exposure 
durations  approximate  thermal  isolation  since  the  heat 
does  not  have  sufficient  time  to  be  conducted  away.  In¬ 
deed,  for  the  cornea,  a  characteristic  conduction  time 
associated  with  CO;  laser  radiation  has  been  given  as 
approximately  10  4  s.“  This  time  is  consistent  with  the 

"M  L-  Wolbarsht,  “Laser  Surgery:  (  O,  or  HF."  ti  ll  ./  Quan¬ 
tum  Electron.  QK-20,  1427-1432  (1984). 
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Figure  27  The  data  points  are  the  experimental  peak 
irradiances  required  to  produce  threshold  damage  as  a 
function  of  beam  diameter.  The  exposure  durations  were 
(a)  0.1  s,  (b)  1.0  s,  and  (c)  10.0  s.  The  curves  are  the  cal¬ 
culated  peak  irradiances  needed  to  produce  a  fixed  tem¬ 
perature  rise,  ATC,  at  a  point  on  the  beam  axis  10  jim 
deep  into  the  epithelium.  ATC  =  41  °C  in  (a),  39°C  in  (b). 
and  36°C  in  (c). 


apparent  asymptotic  limit  suggested  in  Fig.  28.  For  the 
exposure  durations  in  the  present  experiments,  some  of 
the  energy  would  flow  out  because  of  conduction  into 
the  stroma.  Thus  more  energy  would  have  to  be  supplied 
to  compensate  for  energy  lost  through  conduction. 
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Figure  28  Energy  density  at  the  epithelial  damage 
threshold,  ED,h,  for  single-pulse  exposures  as  a  func¬ 
tion  of  exposure  duration.  EDm  appears  to  asymptoti¬ 
cally  approach  a  fixed  value  at  very  short  exposures. 

We  have  incorporated  these  ideas  into  a  new  physi¬ 
cal  model  to  explain  thermal  damage.  Because  of  its 
importance,  we  have  chosen  to  discuss  it  separately  in 
Section  8.0. 

7.2  MULTIPLE-PULSE  EXPOSURES 

The  specification  of  exposure  conditions  for  multiple- 
pulse  exposures  is  considerably  more  complex  than  for 
single  pulses.  In  addition  to  the  peak  irradiance  along 
the  beam  axis,  /0,  the  1/e  diameter  of  the  beam,  d, 
and  the  duration  of  the  individual  pulses,  r  (all  of  which 
were  required  to  specify  the  single-pulse  experiments), 
the  total  number  of  pulses,  N,  and  the  pulse  repetition 
frequency,  PRF,  must  also  be  specified. 


In  the  first  series  of  experiments  that  we  performed 
(Series  1),  we  held  70,  d,  N,  and  PRF  fixed,  while  t 
was  varied  to  determine  the  damage  threshold.  The  pro¬ 
cedure  was  repeated  as  I0,  N,  and  PRF  were  varied 
systematically.  The  results  of  these  experiments  are  listed 
in  Table  2,  where  it  is  apparent  that  threshold  damage 
occurs  when  the  maximum  temperature  increase  is  ap¬ 
proximately  45  °C.  The  temperature-time  histories  that 
correspond  to  the  experimental  parameters  show  a  wide 
variability,  as  indicated  in  Fig.  29.  The  illustrated  tem¬ 
peratures  are  those  on  the  beam  axis  10  /im  below  the 
anterior  surface  and  have  been  calculated  using  an  ex¬ 
tension  of  the  time-dependent  Green’s  function  solu¬ 
tion  discussed  in  Section  3.0. 23  Curves  are  presented 
for  pulse  trains  containing  between  one  and  eight  sub¬ 
threshold  pulses  having  peak  irradiances  of  - 10 
W/cm2.  The  time  histories  for  these  exposures  are 
sufficiently  varied  that  they  can  be  used  to  rule  out  the 
damage  integral  model  (at  least  in  the  form  of  Eq.  1). 
Although  that  model  could  be  used  to  correlate  dam¬ 
age  from  single  pulses,  it  does  not  apply  for  multiple- 
pulse  exposures  because  no  single  set  of  values  of  the 
parameters  A  and  B  can  be  found  for  these  exposure 
conditions. 

Most  of  our  experiments  with  multiple  pulses,  how¬ 
ever,  have  concentrated  on  a  different  experimental  pro¬ 
tocol  (Series  2).  In  Series  2  we  held  d ,  r,  and  PRF  fixed. 
As  the  number  of  pulses  in  the  sequence  was  varied, 
we  determined  the  minimum  value  of  peak  irradiance 
required  to  produce  damage.  Thresholds  were  deter¬ 
mined  in  this  manner  for  several  values  of  PRF  and 
for  two  values  of  r.  The  results  are  compiled  in  Table 
3,  together  with  the  computed  value  of  the  peak  tem¬ 
perature  increase  on  the  beam  axis  10  ^m  beneath  the 
anterior  surface. 


Table  2 

Multiple-pulse  exposures  (Series  1). 


K 


B< 


Figure  29  Calculated  temperature-time  histories  at  a 
depth  of  10  ptm  along  the  beam  axis  for  the  multiple- 
pulse  damage  thresholds  listed  in  Table  2. 
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Table  3 

Multiple-pulse  exposures  (Series  2). 


Plot  Symbol 


PRF 

(Hz) 

Number  of 
Pulses 

7°  2 

(W/cm2) 

Individual  Pulse 
Duration  (s) 

EDlh 

(J/cm2) 

in  Figs. 

30,  31,  and  32 

d 

(cm) 

A  7^ 

CO  ■ 

1* 

1 

86.2 

9.2 

X 

10  3 

0.79 

O 

0.252 

41.4 

1 

4 

79.6 

10.1 

X 

10' 3 

3.21 

O 

0.172 

43.4 

1 

32 

66.8 

9.3 

X 

10'3 

19.9 

o 

0.250 

38.4 

10 

4 

60.0 

9.7 

X 

10  3 

2.32 

□ 

0.256 

43.3 

10 

32 

32.7 

10.9 

X 

I0'3 

11.4 

□ 

0.244 

44.7 

10 

128 

25.3 

9.2 

X 

10-3 

29.8 

□ 

0.270 

39.8 

20 

4 

55.3 

10.0 

X 

10  “3 

2.21 

0 

0.252 

46.4 

20 

32 

27.8 

9.9 

X 

10'3 

8.82 

O 

0.240 

47.7 

20 

128 

16.4 

8.9 

X 

10'3 

18.7 

O 

0.276 

38.8 

1* 

1 

624 

0.96 

X 

10“3 

0.599 

A 

0.102 

54.0 

20 

10 

239 

0.94 

X 

10'3 

2.24 

A 

0.186 

34.7 

20 

100 

171 

0.95 

X 

10-3 

16.3 

A 

0.174 

39.0 

20 

500 

122 

0.92 

X 

10'3 

56.0 

A 

0.190 

33.0 

100 

10 

153 

0.96 

X 

10'3 

1.47 

▲ 

0.202 

36.0 

100 

100 

59.6 

0.95 

X 

10-3 

5.67 

▲ 

0.1% 

35.2 

100  999 

From  Table  1. 

42.0 

0.96 

X 

10'3 

40.3 

A 

0.174 

40.1 

*  Maximum  temperature  rise  calculated  on  the  beam  axis  10  jim  beneath  the  cornea's  anterior  surface. 


Because  of  the  complexity  of  the  temperature-time 
histories  for  experiments  with  multiple  pulses,  it  is  not 
as  straightforward  a  problem  to  deduce  a  general  dam¬ 
age  correlation  as  vs  as  done  for  experiments  with  single¬ 
pulse  exposures.  In  those  experiments,  the  temperature 
increases  monotonically  with  time  until  (depending  on 
the  particular  thermal  model)  damage  occurs  or  the  crit¬ 
ical  temperature  of  the  phase  transition  is  reached  (see 
Section  8.0).  With  multiple-pulse  exposures,  the  tem¬ 
perature  cycles  up  and  down  with  the  pulses,  reaching 
a  somewhat  higher  value  after  each  successive  pulse. 
Still,  it  is  apparent  that  the  maximum  temperature  in¬ 
creases  reached  at  the  end  of  the  pulse  sequence  are  all 
at  least  close  to  the  values  associated  with  damage  in 
the  single-pulse  experiments.  As  noted  above,  this  ob¬ 
servation,  and  the  extreme  variability  of  the  tempera¬ 
ture-time  histories,  rules  out  the  possibility  that  the 
damage  could  be  correlated  by  the  damage  integral  mod¬ 
el  (recall  Eq.  1)  since  it  is  impossible  to  find  a  single 
set  of  parameters  to  satisfy  the  condition. 

Damage  from  multiple-pulse  exposures  to  the  retina 
has  been  examined  at  Letterman  Army  Institute  of  Re¬ 


search,  :j:'  where  it  was  found  that  EDlh  <x  V  .  This 
empirical  relationship  is  examined  in  Fig.  30  for  the  data 
in  Table  3.  The  values  of  ED,h  are  plotted  as  a  func¬ 
tion  of  the  number  of  pulses  on  a  logarithmic  plot.  A 
line  having  slope  ’  i  is  shown  for  comparison.  Although 
there  is  wide  scatter,  the  data  can  be  correlated,  at  least 
roughly,  in  this  manner. 

Since  the  contract  period  covered  by  this  report,  we 
have  reanalyzed  the  data  in  order  to  determine  w  hether 
the  damage  is  correlated  by  an  approximate  critical  tem¬ 
perature  law  or  by  some  form  of  the  modified  critical 
temperature  law.  In  Fig.  31  we  have  plotted  In  Tmdx 
(where  Tmax  =  +  35 °C)  as  a  function  of  the 

total  duration  of  the  pulse  train  for  each  of  the  multiple- 


"4B  F..  Stuck,  I).  .1  Lund,  and  t  S  Beatrice,  Repetitive  Pulse  l  a- 
ser  Data  and  Permissible  Exposure  I  units.  Institute  Report  No 
SR,  1  cltcrman  Arms  Institute  of  Research,  Presidio  ol  San  I  ran 
cisco,  Calif  (Apr  1978). 

*'D.  J  1  und,  B.  P  Stuck,  and  I  S  Beatrice,  Riologtcal  Research 
in  Support  of  Project  MILES,  Inst il ulc  Report  No  %,  I  etter 
man  Army  Institute  ot  Research.  Presidio  of  San  Francisco.  (  aid 
(lul  1981) 
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Figure  30  Energy  density  at  threshold,  EDlh,  as  a 
function  of  the  number  of  pulses,  N,  for  the  data  in  Ta¬ 
ble  3.  The  line  has  slope  V«,  suggesting  that  the  multiple- 
pulse  data  are  approximately  correlated  by  a  relation¬ 
ship  ED,n  <*  A/v‘.  A  similar  correlation  has  been  sug¬ 
gested  for  a  multiple-pulse  damage  to  the  retina.242* 


Exposure  duration  (s) 

Figure  31  This  plot  is  similar  to  Fig.  26  except  that  the 
data  are  for  the  multiple-pulse  exposures  in  Tables  2  and 
3.  Here  the  log  of  the  temperature  on  the  beam  axis  10 
jim  beneath  the  cornea's  anterior  surface  is  plotted  as 
a  function  of  the  total  duration  of  the  pulse  train.  The 
±  10%  error  bars  reflect  both  the  bracket  spacing  in  de¬ 
termining  the  threshold  and  the  estimated  accuracy  of 
the  power  measurement  Unlike  the  single-pulse  data 
in  Fig  26,  there  is  no  evident  trend  with  exposure  dura¬ 
tion  The  data  scatter  about  the  mean  critical  tempera¬ 
ture,  =  75”  (corresponding  to  a  mean  temperature 
increase  of  40  ±  4”C). 


0  1  10  100  100C 

Duration  of  pulse  train  (s) 

Figure  32  Energy  density  at  threshold,  EDm,  for  the 
multiple-pulse  exposures  in  Tables  2  and  3  as  a  func¬ 
tion  of  the  total  duration  of  the  pulse  train.  For  durations 
longer  than  0.1  s,  the  data  are  correlated  by  an  empiri¬ 
cal  power  law  of  the  form  ED,h  ~  D066.  The  plot  sym¬ 
bols  correspond  to  the  respective  conditions  defined  in 
Tables  2  and  3.  The  solid  circles  are  for  the  single-pulse 
exposures  that  have  durations  of  greater  than  0.1  s  (cf. 
Table  1). 


pulse  experiments  in  Tables  2  and  3.  The  total  dura¬ 
tion  of  the  pulse  train  is  given  by  D  =  (N-  \  )/PRF  +  r. 
While  there  is  some  scatter  of  the  values  of  In  Tmix 
about  the  mean,  there  is  no  trend  of  the  data  with  D. 
This  suggests  that,  to  first  order  at  least,  multiple-pulse 
damage  thresholds  can  be  correlated  by  a  critical  temper¬ 
ature  law  in  which  damage  is  associated  with  a  mean 
temperature  increase  of  40  ±  4°C  (7^  =  75  ±  4°C). 
The  correlation  excludes  the  result  for  the  single  1-ms 
pulse. 

A  different,  but  still  empirical,  correlation  of  the 
multiple-pulse  data  is  given  in  Fig.  32.  This  method, 
also  developed  after  the  contract  period  covered  by  the 
present  report,  is  included  here  for  completeness.  For 
durations,  D,  greater  than  ~0. 1  s,  the  energy  density 
at  threshold,  ED,h,  is  seen  to  obey  a  power  law  of  the 
form 


ED,h  =  (const)  D" 


This  correlation  depends  only  on  measured,  not  calcu¬ 
lated,  parameters.  It  is  not  apparent  why  the  points  rep¬ 
resented  by  open  circles  (corresponding  to  the  experi¬ 
ments  for  which  PRE  -  I  Hz,  r  *  10  ms.  and  .V  = 
1,  4,  and  32)  do  not  lie  on  the  same  curse  as  all  the 
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other  data.  However,  the  two  points  for  which  D  > 
0. 1  s  lie  parallel  to  the  straight  line  formed  by  the  rest 
of  the  data,  indicating  that  the  power-law  dependence 
is  the  same.  We  also  have  included  EDth  for  the  single¬ 


pulse  exposures  having  durations  of  >  0. 1  s  on  the  plot 
(solid  circles).  (These  points  are  the  same  as  in  Fig.  28.) 
It  is  evident  that  they  cluster  with  the  multiple-pulse  data 
and  appear  to  be  asymptotic  to  the  D0*6  line. 


8.0  THERMAL  DAMAGE  MODEL 


In  the  section  above,  we  noted  that  certain  aspects 
of  the  epithelial  damage  data  for  single-pulse  exposures 
are  consistent  with  the  hypothesis  that  damage  is  as¬ 
sociated  with  an  endothermic  phase  transition  (see  the 
discussion  pertaining  to  Fig.  28  in  subsection  7.1).  We 
have  incorporated  this  idea  into  a  new  thermal  model 
that  provides  a  physical  basis  for  correlating  epithelial 
damage  data  from  single  pulses.  This  model  explains 
the  slight  variation  in  the  CPT  that  was  noted  by  Eg¬ 
bert  and  Maher4  in  terms  of  heat  conduction  into  the 
cornea  during  the  phase  transition. 

We  consider  the  one-dimensional  thermal  model 
shown  in  Fig.  33,  in  which  a  uniform  flux  of  radiation, 
/„,  is  incident  on  a  layered  cornea  whose  initial  tem¬ 
perature  is  T„.  The  model  is  simplified  further  by  as¬ 
suming  that  the  incident  flux  is  absorbed  in  an  infinites¬ 
imally  thin  section  at  the  anterior  surface  of  the  epitheli¬ 
um  (the  tear  film  is  neglected  in  this  simple  model).  If 
the  thermal  properties  of  the  various  corneal  layers  were 
identical,  then  this  problem  would  be  that  of  the  clas¬ 
sical  semi-infinite  slab  whose  surface  temperature  in¬ 
creases  as  the  square  root  of  time  as  illustrated  in  Fig. 
34a.  The  amount  of  heat  conducted  into  the  inner  layers 
of  the  cornea  depends  on  the  derivative  of  temperature 
with  respect  to  depth,  which,  for  the  simple  slab,  ad¬ 
justs  itself  so  that  all  of  the  incident  flux  is  transported 
into  the  cornea. 

Our  new  model  differs  from  the  simple  slab  by 
postulating  that  an  endothermic  phase  transition  occurs 
once  the  anterior  surface  reaches  the  transition  temper¬ 
ature  denoted  by  T  in  Fig.  34b.  The  time  at  which  the 
front  surface  first  reaches  the  temperature,  Tt ,  is 
denoted  by  I,  .  Just  as  a  water-ice  mixture  stays  at  0°C 
during  melting,  the  surface  temperature  stays  at  7j 
while  the  postulated  endothermic  transition  occurs  (see 
Fig.  34b).  As  a  working  hypothesis,  we  assume  that  ob¬ 
servable  damage  is  associated  with  a  specific  amount 


Epithelium 


Endothelium 


Aqueous 


Phase  transition  in 
the  surface  absorber 

Figure  33  The  geometry  of  a  thermal  model  in  which 
a  uniform  heat  flux  is  incident  on  a  surface  absorber  that 
can  undergo  an  endothermic  phase  transition. 

of  material  per  unit  surface  area  undergoing  the  tran¬ 
sition,  i.c.,  after  a  specific  quantity  of  heat  per  unit  sur¬ 
face  area.  ,  is  absorbed  in  the  phase  transition. 

For  an  exposure  duration,  t  (greater  than  t,  ),  the 
amount  of  heat  per  unit  surface  area,  q,  that  actually 
goes  into  the  transition  is  obtained  by  integrating  the 
difference  between  the  incoming  flux,  and  the 
quantity  of  heat  per  unit  area  that  is  conducted  into 
the  cornea  (  K  dT/dx  |„)  from  time,  t  =  r,  ,  to  r. 
Here  A  is  the  thermal  conductivity  and  the  derivative 


£ 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL.  MARYLAND 


is  evaluated  at  the  front  surface.  For  this  particular  mod-  q  =  j[l  +  ft2] (tan  [6)  -  d J 

el,  the  integral  can  be  done  analytically,  yielding 

x  J[tf  (7;  -r„)]2!/[2F0K]  ,  (5) 


Time  t  (s) 

Figure  34  The  increase  in  surface  temperature  as  a 
function  of  time  for  the  one-dimensional  model  of  a  sur¬ 
face  absorber,  (a)  No  phase  transition  takes  place  and 
the  temperature  increases  as  t (b)  The  same  as  (a)  ex¬ 
cept  that  an  endothermic  phase  transition  with  a  criti¬ 
cal  temperature  of  33*C  above  ambient  occurs.  The 
surface  temperature  would  rise  as  fv’  until  fc  and  then 
stay  constant  until  the  latent  heat  is  supplied. 


where  6  =  [(t -tc)/tc] /!,  and  k  is  the  thermal  diffu- 
sivity  (i.e.,  k  =  K/pc,  where  p  is  the  density  and  c  is 
the  heat  capacity). 

This  model  differs  from  the  one-dimensional  models 
that  have  been  developed  to  examine  laser  damage  in 
solid  state  devices.26.  Those  models,  like  ours,  ac¬ 
count  for  heat  conduction  while  the  temperature  of  the 
material  is  increasing.  However,  they  account  for  the 
endothermic  phase  transition  by  adding  the  energy 
necessary  to  vaporize  the  material  instantaneously;  thus, 
they  ignore  heat  conduction  occurring  during  the  phase 
transformation.  In  applying  our  new  model,  as  discussed 
below,  we  find  that  for  exposures  in  the  1-ms  range, 
the  time  required  for  the  phase  transition  to  occur  can 
be  as  much  as  60%  of  the  total  exposure  duration.  Thus 
the  simplification  that  would  be  afforded  by  assuming 
an  instantaneous  phase  transition  is  not  justified  in  the 
analysis  of  our  experiments. 

If  the  transition  temperature,  Te ,  and  heat  of  tran¬ 
sition,  Qc,  were  known,  then  Eq.  5  could  be  used  to 
predict  the  thermal  flux  of  a  uniform  intensity  beam 
that  would  produce  threshold  damage  for  a  given  pulse 
width,  r.  However,  since  Qc  and  Tc  are  not  known, 
we  have  elected  to  determine  whether  reasonable  values 
for  them  could  lead  to  predicted  values  of  thermal  flux 
and  pulse  width  that  agree  with  the  measured  values 
of  Fn  and  r  at  the  damage  threshold.  Another  diffi¬ 
culty  in  applying  this  model  is  that  the  experiments  are 
conducted  with  a  laser  beam  having  a  Gaussian  inten¬ 
sity  profile;  therefore,  the  damage  would  occur  first  in 
a  region  centered  on  the  beam  axis,  and  the  radial  heat 
conduction  would  remove  some  energy  from  this  re¬ 
gion.  In  the  absence  of  phase  transitions,  radial  heat 
flow  would  cause  the  peak  temperature  to  be  lower  than 
it  would  be  if  the  beam  were  uniform,  as  implied  by 
the  results  on  beam  diameter  dependence  discussed  in 
the  previous  section  (cf.  Fig.  25).  As  a  first  approxi¬ 
mation  to  incorporate  the  effect  of  radial  heat  flow  into 
the  uniform  beam  problem,  we  replace  the  actual  Gauss¬ 
ian  beam  problem,  in  which  the  peak  irradiance  is  /„, 
with  a  uniform  beam  problem  having  an  effective  heat 
flux,  F„.  The  value  of  Fu  is  chosen  such  that,  in  the 
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absence  of  a  phase  transition,  the  calculated  peak  tem¬ 
perature  at  the  end  of  the  exposure  would  be  the  same 
as  it  would  be  for  the  actual  Gaussian  beam.  Having 
done  this,  the  data  points  in  Fig.  35  are  the  experimen¬ 
tal  values  of  the  effective  fluxes  (adapted  from  Table 
I)  that  produce  threshold  damage  at  the  given  pulse 
widths.  The  curve  gives  the  theoretical  predictions  of 
Eq.  5  for  a  transition  that  occurs  at  a  temperature  in¬ 
crease  of  33°C  (T,  =  68°C),  with  Qc  =  0.084  J/cm2. 
If  a  reasonable  fraction  of  the  epithelial  material  in  the 
region  of  the  wound  (e.g.,  5%)  undergoes  the  phase 
change,  then  the  value  of  Qc  would  correspond  to  a 
latent  heat  of  approximately  80  cal/g.  These  values  of 
Tt  and  latent  heat  are  both  physically  reasonable. 
Thus,  in  spite  of  the  simplistic  nature  of  this  new  mod¬ 
el,  these  results  suggest  that  endothermic  phase  transi¬ 
tions  can  explain  the  observed  damage  thresholds. 
Similar  calculations  should  be  performed  for  more 
realistic  phase  transition  models  that  also  account  for 
the  Gaussian  irradiance  profile  and  for  multiple-pulse 
exposures. 


1  300- 


Exposure  duration  |s) 

Figure  35  The  data  points  are  the  experimental  values 
of  the  effective  flux,  F0  (adapted  from  Table  1),  which 
produces  threshold  damage  at  the  indicated  single-pulse 
exposure  duations.  The  ±10%  error  bars  reflect  both 
the  bracket  spacing  in  determining  the  threshold  and 
the  estimated  accuracy  of  the  power  measurement.  The 
curve  is  the  theoretical  prediction  from  Eq.  5  for  the  in¬ 
dicated  parameters.  The  figure  suggests  that  all  of  the 
single-pulse  data  can  be  explained  by  the  model  that 
incorporates  a  phase  transition. 
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